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Preface 


This  independent  study  began  as  an  effort  to  recreate 
and  document  the  most  widely  used  analytical  fallout  code, 
WSEG-10.  Specifically,  local  access  to  this  model  in  computer 
coded  form  will  provide  a  basis  for  further  fallout  studies 
at  the  Air  Force  Institute  of  Technology.  Additional  analyses 
of  the  crossrange  dispersion  term  (o^)  and  model  conservation 
of  activity  were  also  performed.  All  computer  work  was  done 
using  the  ASD  CYBER  74  computer  at  Building  640. 

This  author  gratefully  appreciates  the  guidance  provided 
by  Dr.  C.  J.  Bridgman,  thesis  advisor,  to  accomplish  this 
independent  study.  Thanks  are  also  extended  to  Mr.  Ralph 
Mason,  National  Military  Command  Support  Center,  for  providing 
the  most  recent  computer  coded  version  of  WSEG-10  along  with 
sample  results.  Special  thanks  are  also  extended  to  my  wife 
for  her  patience  through  this  independent  study  and  her  assis¬ 
tance  preparing  this  document. 


Dan  W.  Hanifen 
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Abstract 


The  purpose  of  this  independent  study  is  to  recreate  and 
document  the  most  popular  analytical  fallout  model  in  use  over 
the  past  twenty  years,  WSEG-10.  Local  access  to  WSEG-10  at 
the  Air  Force  Institute  of  Technology,  School  of  Engineering, 
will  provide  a  basis  for  future  fallout  studies.  As  such,  this 
study  provides  a  fully  documented  Fortran  computer  code  contain¬ 
ing  the  most  recent  version  of  the  WSEG-10  analytical  model  with 
sample  output.  To  further  understand  this  computer  code,  a 
general  discussion  of  the  WSEG-10  fallout  model  and  analysis 
of  crossrange  dispersion  (a^)  and  activity  conservation  is 
included.  Results  of  the  analysis  of  a demonstrate  that 
diffusive  growth  is  not  accounted  for  in  the  model  and  that 
crosswind  shear  is  the  dominant,  long  term  effect.  In  a  com¬ 
parative  conservation  analysis,  the  WSEG  model  in  use  today 
does  not  conserve  activity  due  to  the  unnormalized  character 
of  the  crossrange  transport  function.  This  effect  is  sub¬ 
stantial  at  yields  less  than  .1  MT.  Activity  not  conserved 
varied  between  31.4%  at  1  KT  and  a  wind  of  60  st.  mi.  to  less 
than  1%  at  100  MT  and  winds  of  60  st.  mi.  Also  included  is 
a  further  discussion  of  model  limitations  or  inconsistencies 
discovered  either  through  computer  use  during  this  independent 
study  or  during  initial  literature  search. 
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DOCUMENTATION  AND  ANALYSIS  OF  THE 


WSEG-10  FALLOUT  PREDICTION  MODEL 

I .  Introduction 

This  thesis  examines  the  most  popular  analytical  fallout 
prediction  model  in  use  over  the  past  twenty  years  —  WSEG-10. 
The  specific  purpose  of  this  thesis  is  to  recreate  and  docu¬ 
ment  a  working  copy  of  the  latest  computer  coded  version  of 
WSEG-10  (hereafter  referred  to  as  WSEG) .  Local  access  of 
WSEG  will  provide  a  basis  of  comparison  for  future  fallout 
modeling  at  the  Air  Force  Institute  of  Technology  (AFIT) . 

As  such,  the  topics  presented  in  this  thesis  are: 

1.  A  description  of  the  analytical  relationships 
comprising  the  basis  of  the  WSEG  model. 

2.  An  analysis  of  crossrange  dispersion  (o^) . 

3.  An  analysis  of  conservation  of  activity  using  the 
WSEG  model. 

4.  Computer  implementation  using  Fortran  computer 
language  and  the  ASD  CYBER  74  computer. 

5.  A  discussion  of  WSEG  model  limitations  and  short¬ 
comings  . 

The  most  recent  version  of  the  WSEG  model  in  the  form 
of  a  coded  Fortran  subroutine  and  sample  results  were  supplied 
by  Mr.  Ralph  Mason,  National  Military  Command  Support  Center. 
The  subroutine  and  results  are  contained  in  Appendix  A. 
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Since  WSEG  is  based  on  empirical  approximations,  deriva¬ 
tions  of  the  basic  mathematical  relationships  is  often  impos¬ 


sible.  The  following  chronology  is  provided  as  background 
to  the  evolution  of  WSEG. 

Background 

After  the  Mike  nuclear  test  of  the  Ivy  Series  in  October 
of  1952,  the  Rand  Corporation  was  contracted  to  begin  nuclear 
fallout  studies  (Ref.  6:6).  Rand  produced  several  fallout 
prediction  techniques  beginning  with  a  hand-calculated  "Disk- 
tosser"  model  which  was  eventually  converted  to  computer 
solution.  This  early  attempt  was  both  complex  and  time  con¬ 
suming.  It  divided  the  nuclear  cloud  into  many  stacked  disks, 
each  of  which  were  transported  independently.  The  experimental 
data  base  was  poor  and  not  properly  documented.  The  experi¬ 
mental  data  was  available  but  depended  on  such  a  number  of 
parameters  (yield,  height  of  burst,  wind  condition,  shear  and 
soil)  that  generalizations  were  difficult  if  not  impossible. 

As  nuclear  testing  continued,  the  complexity  and  cost  to  re¬ 
fine  and  operate  this  model  grew  along  with  the  data  base. 

As  results  became  reliable,  a  second  technique  was  de¬ 
veloped  by  Rand  to  eliminate  the  costly,  time  consuming,  ma¬ 
chine  calculations.  This  technique  used  analytical  approxi¬ 
mations  to  the  complex  "Disk-tosser"  model.  These,  in  effect, 
considered  the  nuclear  cloud  homogenous  and  "smeared"  the 
fallout  on  the  ground  according  to  the  initia^  parameters  of 
yield,  wind,  shear,  height  of  burst  and  soil  conditions.  The 


resulting  empirical  model  (RM  2460)  used  a  log-normal  activity 
size  distribution  with  a  mean  of  44  microns,  a  standard  de¬ 
viation  of  2  microns,  the  t-i ' 2  (Way-Wigner)  decay  law,  and 
assumed  80%  of  the  activity  is  deposited  locally  (Ref.  13) . 

They  also  calculated  a  function  ^'(T)  representing  the  frac¬ 
tional  rate  of  activity  deposition  everywhere  as  a  function 
of  particle  fall  time  (T) .  Rand  did  not  assign  ^'(T)  a  sin¬ 
gle  functional  type,  although  when  plotted  versus  T,  the 
curve  resembled  a  log-normal  distribution  function  (Refs. 
10:36-40;  4:30).  The  concept  appeared  promising,  but  the 
empirical  model  generated  by  Rand  was  never  popular  (Ref. 

6:13)  . 

In  the  late  1950's,  the  Weapons  System  Evaluation  Group 
(WSEG)  sought  to  create  an  inexpensive,  easy  to  use,  analyti¬ 
cal  fallout  prediction  code  of  their  own.  It  was  published 
in  1959  (Ref.  1).  The  authors,  Pugh  and  Galiano,  incorporated 
Rand  data  for  particle  size  distribution  and  particle  fall 
rates  into  their  original  model  (Ref.  1:27),  They  also 
adopted  or  rediscovered  the  Rand  '  (T)  ,  calling  it  ”g(t)". 

The  WSEG  "g(t)"  represents  the  normalized  fractional  rate  of 
activity  deposition  everywhere  as  a  function  of  time.  It 
was  arbitrarily  assigned  a  negative  exponential  form  which 
empirically  fit  the  Rand  data  everywhere  but  at  very  early 
times  (Ref.  6:13).  In  1960  the  exponent  of  "g(t)"  was  modified 
to  improve  low  yield  capability  (Ref.  2).  In  1962,  a  National 
Academy  of  Sciences  committee  revised  the  WSEG  model  to  its 


present  form  (Ref.  10) .  These  modifications  made  the  model 
more  closely  conform  to  the  experimental  data  collected  during 
the  extensive  nuclear  testing  of  that  decade. 

A  detailed  explanation  of  the  current  form  of  WSEG  is 
contained  in  the  next  section. 
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As  stated  earlier,  WSEG  is  an  empirical  approach  to 
local  fallout  prediction  based  upon  early  nuclear  test  data. 

It  is  designed  to  provide  reliable  fallout  prediction  for 
yields  between  1  KT  and  100  MT.  All  activity  for  deposition 
is  assumed  within  the  fallout  cloud  and  80%  is  assumed  de¬ 
posited  locally.  WSEG  neglects  the  induced  activity  of  the 
stem  created  by  torroidal  circulation  during  the  cloud  rise 
in  early  times  after  the  burst. 

The  version  used  for  this  thesis  assumes  a  land-surface 
burst.  No  adjustments  are  made  to  account  for  burst  heights 
greater  than  zero. 

In  order  to  present  a  description  of  WSEG  in  some  logical 
order,  the  model  will  be  defined  within  a  chronology  of  events 
for  a  nuclear  burst  beginning  with  nuclear  cloud  formation. 

The  cloud  is  initially  formed  because  the  nuclear  fire¬ 
ball  vaporizes  both  the  surface  of  the  earth  at  ground  zero 
and  the  weapon  itself.  The  activity  contained  in  the  cloud 
is  both  neutron  induced  and  fission.  After  formation,  the 
fireball  rises  and  begins  to  cool  at  its  outer  edges  faster 
than  the  center  thereby  creating  the  typical  torroidal  cur¬ 
rents  associated  with  the  nuclear  cloud.  WSEG  arbitrarily 
assumes  that  the  cloud  will  rise  to  a  maximum  center  height 
within  fifteen  minutes  and  then  stabilize.  This  stabilized 
cloud  is  modeled  as  a  right  circular  cylinder  as  in  Figure 


At  this  point,  the  cloud  dimensions  and  altitude  are 
fixed  as  reference  values.  Cloud  center  height  (kilofeet) 
is  given  by  WSEG  as: 


Hc  =  44.  +  6.1  £n(yield)  - 


205  |  Jin  (yield)  +  2 . 42  |  (in(yield)  +  2.42)* 


(1) 


The  radioactivity  in  the  stabilized  cloud  is  assumed 
to  be  normally  distributed  in  both  the  vertical  and  horizon¬ 
tal  directions,  namely 


exp-*5 


p(x,y,h) 


*2  -y.2+ 

f-K  *■ 


<h  -  V 

- - — z - 


(  2  TT  ) 


3/2„2 


(2) 


CJ  O. 

o  h 


Unless  otherwise  stated,  yield  has  units  of  megatons  (or  MT) 


where  x  and  y  are  distances  in  the  downwind  and  crosswind 

directions  respectively  and  h  is  the  height  above  ground  in 

kilofeet.  WSEG  defines  a  and  a.  as 

o  h 


a  (st. mi.)  =  exp (0. 70  +  *n(;(ield>  - 
o  i 

3. 25/(4. 0  +  ( In (yield)  +  5.4) 2 )  (3) 

and  a.  (kilofeet)  =  .18  H  (4) 

n  c 


WSEG  further  defines  the  dimensions  of  the  stabilized  cloud 

where  the  cloud  diameter  is  4a  and  the  vertical  thickness 

o 

is  4a^. 


Independent  of  this  spatial  distribution,  the  radioac¬ 
tivity  is  distributed  on  different  sized  particles  by  some 
activity/size  distribution,  A(r).  As  stated  earlier,  this 
activity/size  distribution  is  based  on  Rand  data  and  defined 
as : 


A  (r) 


exp-J*  ( 


Zn  (m) 
72? 


-  £n(r). 

B _ _ 

Br 


2 


(5) 


where  m  =  (44  microns) 

r  =  (particle  radius  in  microns) 

6  =  .690 

As  time  increases,  the  cloud  will  expand  and  move  hori¬ 
zontally  and  fall  vertically  towards  the  earth  as  fallout 
deposition  occurs.  WSEG  assumes  upward  expansion  is  zero. 
This  horizontal  motion  is  due  to  three  forces  assumed  acting 
on  the  cloud. 
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The  first  force  is  torroidal  circulation  which  at  early 
times  causes  the  fallout  particles  to  be  swept  toward  the 
center  of  the  cloud.  The  resulting  fallout  pattern  is  com¬ 
pressed  around  ground  zero  such  that  the  effective  radius  of 
the  pattern  is  one-half  of  the  actual  radius  of  the  cloud. 
Although  not  explicitly  stated  within  Reference  1,  torroidal 
circulation  is  the  dominant  effect  at  early  time.  This  effect 
lessens  as  torroidal  circulation  decreases  over  time  there¬ 
by  allowing  the  cloud  to  grow  radially.  WSEG  arbitrarily  uses 
three  hours  as  the  cutoff  for  any  torroidal  effect. 

The  second  force  acting  on  the  cloud  is  effective  wind 
(Wind) .  WSEG  used  a  single  effective  wind  vector  over  the 
vertical  extent  of  the  cloud.  The  net  effect  is  to  translate 
the  fallout  cloud  downwind.  Winds  used  to  validate  the  ori¬ 
ginal  model  vary  between  0  and  60  knots. 

The  third  force  is  shear  (S  )  which  WSEG  assumes  is 

c 

constant  over  the  vertical  extent  of  the  cloud.  It  is  the 

change  in  direction  of  the  effective  wind  vector  horizontally 

as  a  function  of  altitude.  Vertical  shear  is  neglected. 

The  shear  used  varies  between  .1  and  .6  knots/kilof eet .  The 

effect  of  the  shear  is  to  expand  the  cloud  and  spread  out 

the  fallout  pattern  as  one  would  open  a  fan. 

The  combined  effects  of  torroidal  circulation,  effective 

wind,  and  shear  are  accounted  for  within  the  terms  describing 

downwind  (o  )  and  crosswind  (o  )  dispersion.  The  function 
x  y 

o  is  affected  by  the  effective  wind  and  defined  as: 
x  1 
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(L  2  +  8a  2) 

ox2(st.mi.)2  =  ao2  L  -2  T  2a  2  (6) 

o  o 

where  L  =  Wind  .  T 

o  c 

and  T  =  Time  Constant  = 

c 

H  H 

1.0573203((|§)HC  -  2.5(g§)2)(l  -  .5  exp-f^f)2)  (8) 

and  aQ  is  defined  by  Equation  (3) .  Wind,  torroidal  growth, 
and  shear  all  affect  ay  which  is  defined  as 

Oy2  (st.mi. )  2  =  aQ2  +  j-  (8|x  +  2ax|ao2)  + 

b  ‘WhV 2  +  b  <  <x  +  2ox>  Vc-W 2  <9) 

where  oq  and  ah  are  defined  by  Equations  (3)  and  (4)  re¬ 
spectively,  Sc  is  the  crosswind  shear  and 

L2  =  L  2  +  2a  2  (10) 

O  X 

Thus,  ax  is  fixed  by  initial  cloud  parameters  and  effec¬ 
tive  wind  while  a will  vary  with  time.  This  is  discussed 
in  depth  in  the  next  chapter. 

Throughout  the  growth  and  transport  of  the  radioactive 
cloud  there  is  a  continual  fall  of  particles  back  to  the 
ground.  As  mentioned  in  the  Background,  WSEG  states  that 
there  must  be  some  function  "g(t)"  which  describes  the  frac¬ 
tional  rate  of  activity  arrival  on  the  ground  everywhere  at 
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some  time  t.  The  integral  of  this  function,  G(t),  represents 
the  fraction  of  activity  down  at  time  t  where 

G(t)  =  o/fcg  (f)  dt'  (11) 

This  g(t)  function  will  be  independent  of  the  horizontal 
activity  distribution  and  therefore  independent  of  the  growth 
of  0^  with  time.  On  the  other  hand  g(t)  will  be  dependent 
on  the  initial  vertical  distribution  and  the  activity/size 
distribution  which  determines  particle  fall  rate. 

This  activity  deposition, g (t) ,  is  assumed  in  the  original 
WSEG  document  without  derivation  as 


g(t) 


n  .  t .  n 

P  exp-  (^-)  o 

c 

tc  r(i  +  1  ) 


(12) 


where  T  =  Time  Constant 

c 

nQ  =  1.5  -  ,25("|)!  (13) 

F  =1.0 

This  arbitrary  choice  of  g(t)  is  based  on  Rand  calcula¬ 
tions  which  assume  an  activity/size  distribution  given  by 
Equation  (5).  These  calculations  are  neither  shown  nor  re¬ 
ferenced  in  the  original  WSEG  model  contained  in  Reference 
1.  If  the  activity/size  distribution  for  a  given  set  of 
initial  conditions  is  different  than  that  given  by  Equation 
(5),  the  form  of  g(t)  should  change.  This  is  not  possible 
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under  the  WSEG  model  where  the  function  g(t)  is  fixed  as 
Equation  (13) .  The  only  possible  compensation  for  various 
activity/size  distributions  results  because  Tc  varies  with 
yield  (Ref.  2) . 

After  modifications  in  1962, 


F  exp  -  (^-)n 
c 

t  r  (i  +  l) 

^  n 


where 


n  =  F  =  1.0 
o 


n  L  2  +  a  2 
o  o  x 

L  z  +  .5 a  7 


The  function  g(t)  can  be  transformed  to  a  downwind  distance 


function  since 


g(t)dt  =  g(x)dx 


x  =  Wind  .  t 


The  downwind  distance,  x,  is  substituted  into  Equation  (14) 
for  t  as: 


g(x)  = 


,  t  .  Wind  .  n 

-  wind1 
c 

Tc  .  Wind  .  T(1  +  1) 

n 


g(x)  = 


exp-f^)" 

Lo  "r(i  +  1  ) 
n 
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To  provide  continuity  in  a  "0"  wind  environment  near  ground 
zero  WSEG  replaces  Lq  with  L  and  the  domain  is  arbitrarily 
extended  by  setting  x  =  |x|.  Therefore  the  final  activity 
deposition  function  in  terms  of  distance  is: 


g(x) 


exp-(M,n 

l  r (i  +  l) 


(18) 


n 

g(x)  in  this  form  also  represents  the  fallout  deposition 
distribution  function  used  within  WSEG. 

In  order  to  predict  upwind  fallout  and  at  the  same  time 
preserve  normalization,  a  function  <j>  is  empirically  inserted 
where 


4>  = 


1 

/2tt 


exp(-bz2)dz 


(19) 


and 


x 

o  a 
x 


) 


The  normalized  downwind  and  upwind  distribution  is  then 
represented  as : 


x 

V*1 


)  *g(x) 


(20) 


where 


x 

axai 


)  .  g (x) dx 


1 


The  model  forces  <j>  to  behave  as  follows: 

.5  <  i)  <  1.  for  x  >  0. 
0.  <  *  <  .5  for  x  <  0 . 
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aa  is  an  adjustment  factor  to  reduce  the  area  covered  by  fall¬ 
out  prior  to  cloud  stabilization  due  to  torroidal  compression 
(Ref.  4:35)  and  is  defined  as: 


ota 


_  1 _ 

(1  +  .001  .  Hc  .  Wind) 


(21) 


where  oq  is  defined  by  Equation  (3) . 

The  parameter  "n"  defined  by  Equation  (15)  and  used  as 
the  exponent  in  g(x)  and  get),  is  plotted  in  Figure  2  for 
yields  of  1,  10,  100,  1000  and  10,000  KT.  As  seen,  "n"  is 
a  weak  function  of  yield  and  varies  dramatically  for  winds 
between  0  and  1  However,  for  all  intents  and  pur¬ 

poses,  the  variation  can  be  described  as: 

n  =  2  when  Wind  -  0 
n  =  1  when  Wind  >  0 

Figures  3  and  4  depict  g(x)  and  cj>.g(x)  for  a  1  MT  burst. 

Figure  3  represents  the  case  where  the  effective  wind  and 

shear  are  both  0  while  Figure  4  represents  the  behavior  of 

g(x)  and  ij).g(x)  where  the  effective  wind  is  10  and 

shear  is  0.  If  plotted  further  downwind  g(x)  and  4>.g(x) 

would  asymptotically  approach  0.  The  validity  of  the  sub- 
00  00 

stitution  /  <b.g(x)dx  =  1  for  /  g(x)dx  =  1  is  addressed 

—  oo  *  •'  O 

in  Section  IV  concerning  conservation  of  activity.  Numerical 
integration  was  used  to  integrate  g(x)  and  <f>.g(x)  and  as 
such  finite  integration  limits  were  established  as  an  approx¬ 
imation. 
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Figure  2.  Parameter  "n”  vs.  Wind 


'  The  downwind  transport  function  £  can  now  be  written 

as : 

f  =  Yield  .  SNC  .  <J)  .  gCx)  .  fission  fraction  C22) 

where  SNC  is  the  Source  Normalization  Constant  = 

2  x  106  Roentgens/hr/MT/(st.mi. ) 2 

(  The  crosswind  transport  function  is  a  modified  Gaussian  dis¬ 

tribution  of  the  form: 
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(23) 


exp 


-hi-*—) 

‘  V/V  -  ' 


f  = 

y 


a2a 


vmr 


where  a2  is  an  adjustment  factor  added  in  1962  similar  to 
aj  but  only  effective  for  2  hours  (Ref.  4:  36).  it  is 
defined  as: 


a2 


(1  +  .001  .  Hc  .  Wind 


(1  -  4>  ( 


2x 

Wind 


)) 


(24) 


The  activity  of  the  fallout  in  the  cloud  decays  by  Way- 

_i.  2 

Wigner  (t  )  as  already  mentioned,  as  does  the  fallout  de¬ 
posited.  This  assumes  no  fractionation.  WSEG  creates  iso¬ 
dose  rate  contours  by  utilizing  the  Unit  Time  Reference  Dose 
Rate  (dh+^)  which  is  the  product  of  the  downwind  and  cross- 
wind  transport  functions  (f^  .  f^.)  .  represents  the  ac¬ 

tivity  at  some  point  (x,y)  one  hour  after  detonation.  This 
includes  all  activity  that  has  arrived  at  (x,y)  in  1  hour 
plus  all  activity  that  will  be  deposited.  These  contours 
are  elliptical  with  the  major  axis  along  the  hotline.  The 
length  of  the  contour  depends  on  the  initial  yield  of  the 
weapon  and  on  the  magnitude  of  the  effective  wind  vector. 
Contour  width  is  determined  primarily  by  shear  (see  Section 


III)  . 

To  obtain  a  measure  of  dose  to  humans,  "Biological  Dose" 
was  defined  as  the  product  of  the  D^+^  and  a  conversion  fac¬ 
tor,  called  Bio.  Bio  is  an  empirical  function  depending  on 
fallout  arrival  time  and  length  of  exposure.  Ten  percent  of 
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the  dose  received  is  assumed  irreparable  and  ninety  percent 
is  assumed  reparable  with  a  thirty  day  time  constant  (Ref. 
11).  Mathematically  this  set  of  conditions  is  written  as: 

D ( t)  =  10%  t  /tDH+1T-1’2  d-r  + 

90%  t  /tDH+1T  '' 2exp(|(t-T) )dT  (25) 

cl 

where  K  =  30  days  and  t  =  average  expected  time  of  arrival 
of  the  fallout  and  is  defined  as: 


t 

a 


(0.25  + 


L  2 .x  +  2a  2  T 
o  (  x  ) 

L2(L  2  +  .5a  * 
o  x 


c: 

) 


2 


+ 


2a  -  2Ti 2  , 

X  J  .  *5 

L~2+.5a  ' 
o  x 


(26) 


This  equation  assumes  the  earliest  arrival  time  of  fallout 

anywhere  is  .5  hours.  Ti  equals  one  hour  and  is  included 

to  maintain  dimensionality.  It  is  often  eliminated  from  the 

expression.  At  large  x  Equation  (26)  reduces  nicely  to 

t  =  x- 
a  Wind 

Equation  (25)  was  solved  numerically  and  plotted  as 
Dose  vs.  Time.  Bio  was  then  approximated  in  Reference  1  as: 

Bio  =  (t/19)  ” 


so  that  the  dose  at  some  time  after  activity  arrival  is 
defined  as: 


Dose 


H+l 


Bio 


18 


Further  refinements  in  the  model  resulted  in  a  second 
order  approximation  for  Bio  of  the  form: 

t 

Bio  =  exp- (.287  +  .  52  £n  + 

t 

.04475  £n(II^g)2)  (27) 

which  is  in  use  today. 

These  special  conditions  dictated  the  necessity  for  the 
definition  of  a  special  unit  of  dose,  the  ERD.  The  ERD  or 
Equivalent  Residual  Dose,  actually  has  units  of  Roentgens 
even  though  it  pertains  to  human  whole-body  damage.  This  is 
not  in  keeping  with  the  current  philosophy  in  assigning  units 
of  exposure  and  dose.  This  use  of  Roentgens  as  a  measure  of 
dose  instead  of  exposure  tends  to  confuse  those  new  to  WSEG. 
The  subroutine  contained  in  Appendix  A  refers  to  the  original 
Pugh  definition  of  dose.  The  AFIT  version  contained  in 
Appendix  B  generates  only  dose  rate  contours.  Conversion  to 
the  proper  units  will  be  necessary  if  dose  contours  are  re¬ 
quired  in  the  future. 


III.  Crossrange  Dispersion 


This  section  will  examine  the  crossrange  dispersion 
( )  of  the  fallout  cloud  by  developing  the  terms  in  its  de¬ 
finition.  The  specific  purpose  of  this  analysis  is  to  deter¬ 
mine  whether  wind  shear  or  torroidal  growth  is  the  most  signi¬ 
ficant  contributor  to  crossrange  dispersion.  Also  an  ex¬ 
amination  of  the  resemblance  between  the  form  of  the  torroi¬ 


dal  growth  term  and  diffusive  growth  according  to  Fick's  Law 
is  included. 


Recall  o  2 

y 


is  defined  by  Equation  (9)  as: 


2  =  o  2  (1  +  1 x  +  2a  I  )  +  (a  2T  2a,S  2)  + 

7  o  L1  x1  L^  xchc 


fr((x  +  2a  )  L  ToS) 2 
L  x  o  c  h  c 


This  discussion  will  first  consider  the  contribution 
to  crossrange  dispersion  by  shear  represented  by  the  second 
and  third  terms  of  Equation  (9) .  The  first  term  expressing 
torroidal  growth  is  discussed  later  in  this  section. 

To  begin  with,  the  second  term  of  Equation  (9)  was  added 
to  the  WSEG  model  without  derivation  or  reference  as  an  after¬ 


thought  to  reduce  the  fallout  concentration  near  ground  zero 
using  shear  (Ref.  10) .  Little  more  can  be  said  about  its 
development.  In  fact,  its  total  contribution  to  a^  is  small 
considering  the  remaining  shear  term. 

This  remaining  shear  term  represents  the  original  shear 


contribution  to 


that  Pugh  and  Galiano  postulated,  modified 
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without  reference  or  derivation  for  better  response  near 
ground  zero  for  low  winds  (Ref.  1  :14).  The  term  can 

easily  be  transformed  to  the  original  form  in  Reference  1 
which  is  (Sco^  ^Ind^2'  T^e  derivation  will  utilize 

this  early  form  and  will  consider  shear  as  the  sole  contrib¬ 
utor  to  crossrange  dispersion,  or  *  Thi-S 

discussion  is  taken  from  pp.  10-35  of  the  original  WSEG 
document  to  derive  the  expression  (Ref.  1) .  WSEG  assumes 
that  effective  wind  is  the  average  of  all  wind  vectors  through 
which  the  fallout  particles  travel  to  earth. 

This  effective  fallout  wind  is  defined  as: 

wind  -  h/0  ^TkWc/  vWdh  (28) 

where  W(h)  is  local  wind  at  altitude  h  and  V(h)  is  the  rate 
of  fall  of  a  typical  particle  at  altitude  hQ.  Since  the 
wind  data  is  obtained  at  discrete  altitudes  it  can  also  be 
represented  as: 

I  T±W(i) 

Wind  =  -  (29) 

n 


where  =  time  spent  within  each  wind  layer  by  a  typical 

•  th 

particle  falling  to  earth  and  W(i)  is  the  wind  in  the  l  layer. 

WSEG  further  assumes  that  the  effective  wind  is  a  slowly 
varying  vector  dependent  upon  altitude  which  can  be  expanded 
around  cloud  center  height  as  a  Taylor  series  for  both 
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crosswind  (W^)  and  downwind  (.W^)  components: 

dW  d2W 

Mx  -  \h  +  (dir>h  (h  -  V  +  *  <h  -  ho»’  +  ••• 

no  o 

o 


My  -  0  +  <dh*>  <h  -  ho>  +  ’‘<diT*h  <h  -  V'  +  •• 


un  , 

with  Downwind  Shear  =  S  =  .  (r — r  •  y11^ — r) 

x  dh  h  hr-kilofeet 

o 


Crosswind  Shear  =  S  =  (-5^-),  (t—  ,  . ,  ~c — r) 

c  dh  h  hr-kilofeet 
o 


The  higher  order  terms  are  neglected  as  they  are  assumed 

small.  WSEG  also  assumes  downwind  shear  is  neglected. 

Finally  Sc  .  =  magnitude  of  wind  vector  direction 

change  over  1  standard  deviation  in  altitude  and  the  total 

contribution  to  a  is  S  a,  t,  where  "t"  is  time  after  burst 

y  c  h 

in  hours  and  defined  as  .  Pugh  and  Galiano  recommend 

using  shear  evaluated  for  layers  to  two  to  four  o.  above 


and  below  H  for  best  results  (Ref.  1:24), 
c 


Figure  5 


demonstrates  how  the  shear  affects  assuming  the  cloud  is 
modeled  as  a  cylinder  with  a  radius  equal  to  the  radius  of 
the  cloud. 

Figures  6-8  show  that  wind  shear  contributes  signifi¬ 
cantly  at  later  times.  These  figures  depict  o^2  versus 
time  for  several  shear  and  wind  conditions.  The  two  terms 
representing  a  shear  contribution  were  combined  into  one 
curve  and  included  on  the  same  graph.  The  torroidal  growth 
term  was  also  included  as  a  separate  curve.  All  curves 
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Yield  =  10  MT 

K.OO’i  Fission  Fraction 

1  Wind  »  20  st. mi . 
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were  generated  using  the  subroutine  Dose  with  an  appropriate 
main  program  designed  for  plotting  with  the  Calcomp  plotter 
in  Building  640.  Additional  graphs  are  available  for  sever¬ 
al  wind  and  shear  combinations  in  Appendix  C. 


Torroidal  Growth 

As  seen  from  Figures  6-8,  the  torroidal  growth  term 
predominates  at  early  time.  It  is  defined  in  Equation  (9) 

Q 

as  o  2(1  +  rU  +  2o  I).  In  the  absence  of  shear,  WSEG 

assumes  this  expression  relates  the  compressed  dimensions  of 

the  fallout  pattern  near  ground  zero  to  the  fallout  cloud 

dimensions  as  a  function  of  time.  WSEG  defines  the  radius 

of  the  fallout  pattern  as  a  ,  the  "effective  radius"  where 

the  effective  radius  is  arbitraily  assumed  equal  to  one-half 

the  radius  of  the  stabilized  fallout  cloud  at  15  minutes  and 

where  o  2  =  a  2  (1  +  |^  |x  +  2a  I). 
e  o  L  1  x 1 

Initially  Pugh  and  Galiano  defined  og2  as  oq2 (1  + 

•  «  X  » 

without  derivation  or  reference  where  represents  time 


after  burst  in  hours  (Ref.  1:13). 


It  is  this  form  which 


resembles  simple  diffusion  according  to  Pick's  Law.  There 
are,  however,  several  serious  inconsistencies  with  this  as¬ 
sertion  which  are  evident  in  the  following  discussion  which 
develops  both  the  present  form  of  ag2  used  in  Equation  (9) 
and  the  diffusivity  parameter  (Dv)  for  Fick's  Diffusive  Law 
(see  Appendix  D) . 

First,  while  the  expression  for  ag2  resembles  diffusive 
growth,  it  was  originally  intended  to  allow  for  torroidal 


27 


IP"*'  1  > 


growth  (Ref.  9) .  In  either  case,  without  modification,  the 

initial  expression  for  a e2  is  dimensionally  incorrect! 

Pugh  and  Galiano  corrected  this  problem  by  substituting 

the  ratio  L/T  for  Wind  where  L  is  defined  as  Wind  .  T  and 
c  o  c 

L  =  L  when  L  2  >>  2a  2  (see  Equations  (7)  and  (10)).  At 
o  o  o 

this  point,  Pugh  and  Galiano  assigned  Tc  a  dimensionless 
value  of  8  which  corrected  the  units  problem  and  compensated 
for  various  yields  greater  than  1  megaton.  There  is  no  ap¬ 
parent  reason  for  deleting  the  units  except  convenience. 

Two  further  modifications  were  made  to  o e2  resulting 
in  its  present  form.  The  first  modification  set  x  =  |x|  to 
account  for  both  downwind  and  upwind  fallout  pattern  growth. 
The  second  modification  set  |x|  =  |x  +  2ox|  which  prevented 
a  minimum  at  ground  zero. 

To  derive  an  expression  for  Diffusivity  (Dv)  according 
to  Fick's  Law,  og2  is  first  defined  as  aQ2(l  +  -^x L)  which 
is  a  good  assumption  if  X  >>  2ax-  Wind  .  Tc  is  substituted 
for  L  by  the  same  reasoning  used  earlier  to  develop  this  ex¬ 
pression  and  the  result  is  oe2  =  •  °02  +  aQ2  = 

c 

O 

a  2  +  o  2t.  T  and  o  are  allowed  to  vary  according  to 

O  T  O  C  O 

c 

yield,  t  is  time  in  hours  after  burst.  Diffusivity  according 
to  Fick's  Law  is  therefore: 

D  =  D  =  ~-o  2  ( (Su •"■-—)  (See  Appendix  D)  (30) 

WSEG  v  T  o  '  hours 
c 
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Assuming  DWSEG  is  dimensionally  correct,  a  comparison 
of  the  Diffusivity  parameter  and  diffusive  growth  was  made 
between  WSEG  and  the  Department  of  Defense  Land  Fallout  Pre¬ 
diction  System  (DELFIC) .  The  following  data  was  supplied  by 
Major  Scott  Bigelow,  Air  Force  Weapons  Lab  (AFWL)  using 
DELFIC  for  a  single  particle  size  group: 

Yield  =  .1  megatons 

H  =9.0  kilometers 
c 

o  =  2607  meters 
o 

Oj  =  2863  meters 

t  =  611  sec 
o 

tf  =  2.6  x  104  sec 

Radius  of  particle  group  =  33.8  microns 

Diffusivity  was  calculated  from  the  above  data  using  DELFIC 
as : 

D  =  5  x  10~ 6  meters* 
v  sec 

where  the  quantity  (c^  -  aQ)  represents  cloud  growth  due  to 
the  diffusive  process  in  time  (t^  -  tQ) . 

For  .1  MT  burst  using  WSEG:  Mean  radius  of  particle  for 
WSEG  activity/size  distribution  =  44  microns. 

T  =  4.995  hours 
c 

Hc  =  9.1  kilometers  (Equation  (1)) 

n  =  4  j J°1  (st  .mi  .  )  2 

WSEG  4.995  hr 
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and 


oq  =  .736  statute  miles  (from  Equation  (3)) 


Therefore 

DWSEG 


4(.736  st. mi.)2  _  4(1182.86  meters) 2 

(4.995)  (3600  sec)  (4.995)  (3600  sec) 
hr  hr 


dwseg  =  310-96 


meters* 

sec 


Clearly  the  magnitude  of  "dif fusivity"  represented  in 
WSEG  is  something  entirely  different  than  the  slow  diffusive 
process  represented  by  the  DELFIC  Dv-  Further,  assuming  the 
fallout  cloud  radius  in  DELFIC  can  be  represented  as  2oe 
(as  in  WSEG) ,  the  radius  shows  a  diffusive  growth  of  512 
meters  in  7.0525  hours.  Using  the  identical  time,  the  WSEG 
cloud  radius  would  show  a  growth  of  1.58  x  10 7  meters.  Ap¬ 
plying  the  cutoff  of  three  hours,  the  growth  is  still  7.72 
x  106  meters! 

Therefore,  while  the  form,  excluding  dimensionality, 
appears  to  fortuitously  resemble  diffusive  growth,  it  does 
not  provide  proper  parameters  as  diffusion  is  a  small  part 
of  cloud  growth.  Also  the  model  places  a  three  hour  time 
limit  on  the  effects  of  this  term  which  corresponds  to  a 
contribution  to  o^2  of  13.334  (st.mi.)2.  This  also  is  not 
reasonable  for  diffusive  growth  as  it  would  continue  until 
fallout  deposition  is  completed. 
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IV.  Conservation 


The  purpose  of  this  section  is  to  examine  the  capability 
of  WSEG  to  deposit  within  the  fallout  pattern  all  activity 
assumed  initially  present  in  the  fallout  cloud  at  t  =  15  min¬ 
utes.  This  total  activity  is  a  product  of  yield.  Source 
Normalization  Constant  (SNC) ,  and  fission  fraction  which  is 
taken  in  WSEG  as  2  x  10 6  (yield) (fission  fraction) 

Roentgens (st.mi. ) 2 
hr 

In  order  to  recover  this  product,  DH+^  was  integrated 
over  the  entire  fallout  pattern.  Crosswind  integration  was 
accomplished  analytically  while  the  downwind/upwind  integra¬ 
tion  was  accomplished  numerically  using  a  trapezoidal  tech¬ 
nique.  Also  included  is  a  discussion  of  the  effect  on  con- 
servation  of  activity  in  WSEG  by  substituting  4>.g(x)dx 
for  Q/ag(x)dx  within  Equation  (21)  where  a  and  b  represent 
finite  integration  limits  used  to  make  the  trapezoidal  in¬ 
tegration  possible. 

This  examination  specifically  focuses  on  the  crosswind 
transport  function,  f^,  defined  by  Equation  (22)  and  its 
effect  on  conservation  as  yield  and  Wind  are  varied.  The 
function  f  is  not  properly  normalized  due  to  the  addition 
of  a2.  In  order  to  evaluate  the  effect  of  f  on  conserva- 

y 

tion  and  eliminate  uncertainty  in  the  results  due  to  the 
numerical  integration  technique  used,  the  recovered  product 
was  compared  to  results  obtained  under  identical  conditions 
in  an  identical  manner  using  the  crosswind  transport  function 
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This 


originally  defined  by  Pugh  and  Galiano  without  a2 
original  transport  function  is  defined  as: 


exp- (^) 


f  = 

y 


(Ref.  1:10) 


For  the  purposes  of  this  section  the  original  transport  func¬ 
tion  will  be  distinguished  from  Equation  (19)  by  referring  to 
it  as  f 

yo 

The  following  subsections  discuss  the  specific  method 
by  which  conservation  was  examined  and  the  results  of  this 
examination.  These  results  are  tabulated  in  Tables  I  through 
III  in  this  section.  All  calculations  were  done  using  the 
ASD  CYBER  74  computer  using  subroutine  Dose. 

Method  /Results 

The  total  activity  within  the  fallout  ground  pattern  is: 


/  Iff  dxdy  =  SNC.  fission  fraction  .  Yield 
— oo  —  oo  x  y  "* 


where  f  and  f  are  defined  by  Equations  (21)  and  (22)  and 
x  y 

fission  fraction  =  1.0.  The  product  f^  .  f  at  any  x  and  y 
defines  DH+^  at  that  location.  The  upwind/downwind  integra¬ 
tion  limits  were  replaced  by  finite  values,  a  and  b,  which 
are  defined  on  page  34 . 

To  reduce  the  numerical  integration  in  Equation  (32)  from 
two  dimensions  to  one,  the  crosswind  integration  was  first 
accomplished  analytically  from  -®°  to  +°°  using  the  properties 
of  a  standard  Gaussian  distribution.  In  this  case  however. 
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CO 

_mf  fydy  /  1  since  f  is  not  properly  normalized  because  of 
02-  Recall  from  Equation  (23)  that  a2  is  a  function  of 
yield,  wind  and  downwind/upwind  distance.  This  dependence 
on  x  is  not  subscripted. 

It  can  be  shown  at  any  x: 


f  -  <*2 

y  a  2 


a,  exp-h  (—-%-) 
f  =  *  V 

y  /2no  c*2 


where 


exp-%  (-  ^----) 
CL2  oy 

v'2-no  a.2 


is  a  standard  normalized  Gaussian  distribution  which  if  in¬ 
tegrated  from  -°°  to  +“  would  equal  1.  Therefore  at  any  x: 

exp-Jj  ( — ^— )  2 
*  a  20 

°o  oo  y 

/  f  dy  =  a  2  /  - Tjt - -  dy  =  1  .  a2 

-oo  y  /2-no  a 2  1 

1  y 


If  y  =  0,a  simple  substitute  for  the  crosswind  integra¬ 
tion  is  available  because  f  (0)  =  -77^ —  and  therefore 

y  /2no 

1  y 

a2fy (0)  /TnOy  =  1  .  az.  Placing  the  relationship  in  Equation 
(27)  reduces  it  to  the  following: 


/a a2f  (0)  a  /2irf  dx  =  SNC  .  Yield 
by  y  x 


(33) 


where  fy(0) 


f  =  dtii,  along  the  hotline.  This  hotline  is 
x  H+l  r 


an  imaginary  line  extending  directly  downwind  or  upwind  of 
ground  zero  where  maximum  activity  is  deposited.  (See  x 
axis  in  Figure  1) .  The  function  o and  02  are  defined  by 
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Equations  (9)  and  (23)  respectively  and  a  and  b  define  the 
finite  integration  limits  explained  in  the  next  paragraph. 


The  shear  contribution  to  a 

y 


was  neglected  throughout  this 


section  as  it  was  found  to  have  no  effect  on  conservation. 


The  remaining  upwind  and  downwind  integration  was  accom¬ 
plished  simply  by  a  trapezoidal  integration  of  DH+^  along 
the  hotline  beginning  at  ground  zero.  Criteria  for  the  in¬ 
tegration  limits  was  based  upon  dose  rate  contours.  Trial 
and  error  determined  the  limiting  contour  that  maximized  the 
recovered  activity  in  a  zero  Wind  condition  where  012  is  in¬ 
effective  while  minimizing  computer  time.  Activity  lost  was 
less  than  .1%  in  the  cases  examined  when  the  .1  Roentgen/hour 
dose  rate  contour  was  used.  Maximum  upwind  (a)  and  downwind 
(b)  distance  traveled  to  the  .1  Roentgens/hr  dose  rate  con¬ 
tour  were  used  as  integration  limits  and  noted  for  comparison. 
These  limits  vary  with  yield  and  wind.  As  a  second  check, 
<f>.g(x)  was  also  trapezoidally  integrated  between  a  and  b  to 
be  sure  that  all  activity  available  for  deposition  was  de¬ 
posited.  Step  size  for  g(x)  and  dose  rate  integration  was 
identical  for  each  wind  and  yield  condition.  These  step  size 
varied  from  .0001  to  .1  statute  miles  depending  on  the  con¬ 
ditions.  The  results  are  tabulated  in  Table  I. 


It  can  be  seen  that  the  present  WSEG  model  is  not  con¬ 
servative  by  some  average  effective  02  where  this  average 
ct2  is  the  ratio  of  the  recovered  activity  to  the  initial 
activity.  Table  I  indicates  a  significant  reduction  in  re¬ 
covered  product  (SNC  .  Yield)  at  low  yields  and  high  Wind 
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TABLE  I 


Comparison  of  Initial  vs.  Recovered  SNC  .  Yield 
For  the  Present  Version  of  WSEG 


Yield 

Wind 

st.mi. 

Initial 

SNC  .  Yield 

R-  { st .mi . ) * 

Recovered 

SNC  .  Yield 

R- (st .mi . ) * 

Cumulative 
g  (X) 

(per  st.mi.) 

hr 

hr 

hr 

1KT 

0 

2  x  103 

2.000 

X 

103 

1.0000 

30 

1.590 

X 

103 

0.9991 

60 

1.370 

X 

103 

0.9983 

10KT 

0 

2  x  10“ 

2.000 

X 

10“ 

1.0000 

30 

1.782 

X 

10" 

0.9997 

60 

1.678 

X 

10" 

0.9994 

100KT 

0 

2  x  10s 

2.000 

X 

10s 

1.0000 

30 

1.931 

X 

10s 

1.0000 

60 

1.889 

X 

10s 

1.0000 

1MT 

0 

2  x  106 

2.000 

X 

10s 

1.0000 

30 

1.970 

X 

106 

1.0000 

60 

1.948 

X 

10® 

1.0000 

10MT 

0 

2  x  107 

2.000 

X 

107 

1.0000 

30 

1.985 

X 

107 

1.0000 

60 

1.973 

X 

107 

1.0000 

100MT 

0 

2  x  108 

2.000 

X 

10® 

1.0000 

30 

1.992 

X 

10® 

1.0000 

60 

1.986 

X 

10® 

1.0000 

*"R"  is  an  abbreviation  for  Roentgens. 
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since  a  2  ^  1.  When  the  Wind  =  0,  the  recovered  activity  at 
all  yields  is  identical  to  the  activity  recovered  using  f 

(defined  by  Equation  (31))  in  a  zero  Wind  condition.  In  this 
condition  a2  =  1.0  and  has  no  effect  thereby  normalizing  the 
crossrange  distribution  function.  The  lower  yields  are  more 
dramatically  affected  since  much  more  activity  is  deposited 
in  a  short  time  period.  Recall  a2  is  effective  for  2  hours 
after  stabilization.  The  net  result  is  to  reduce  the  total 
activity  deposited  within  the  .  1  R— ^g-en  -  dose  rate  contour 
by  the  percentages  shown  in  Table  II.  Cumulative  <{>.g(x)  is 
very  close  to  1.0  in  all  cases  indicating  that  all  activity 
has  been  deposited. 

The  case  of  the  original  distribution,  f^Q  defined  by 
Equation  (31) ,  was  handled  in  the  same  manner  as  described 
earlier  where: 

/°°f  ( 0)  dy  =  1 

-°0  yo  2 

Again,  a  simple  substitute  is  available  to  reduce  the  nec¬ 
essary  integration  to  one  dimension  through  the  properties 

of  a  Gaussian  distribution.  If  y  =  0  then  f  (0)  =  -pr—  -. 

yo  r  2tt o 

Thus  o  /2tt  f  (0)  =  1  and  this  expression  is  substituted  in- 
y  yo 

to  Equation  (32)  which  simplifies  to: 


,Sa f  (0)  a  /2tt  f  dx  =  SNC  .  Yield 
b  y  y  x 


where  f  (0)  .  0  =  D„.,  along  the  hotline  and  0  is  defined 

y  x  H+l  y 

by  Equation  (9) .  The  integration  limits  are  a  and  b. 
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TABLE  II 


Per  Cent  Decrease  of  Recovered  Activity 
vs.  Yield  and  Wind 


Yield 

Wind 
(st.mi. ) 
hr 

Per  Cent 
Decrease 

1KT 

0 

0.0 

30 

20.6 

60 

31.5 

10KT 

0 

0.0 

30 

10.9 

60 

16.1 

100KT 

0 

0.0 

30 

3.4 

60 

5.6 

1MT 

0 

0.0 

30 

1.5 

60 

2.6 

10MT 

0 

0.0 

30 

0.8 

60 

1.4 

100MT 

0 

0.0 

30 

0.4 

60 

0.7 

Again,  upwind  and  downwind  integration  was  accomplished 
by  trapezoidally  integrating  DH+^  along  the  hotline  beginning 
at  ground  zero.  Integration  limits  and  step  sizes  were  iden¬ 
tical  to  those  used  for  the  present  distribution  for  each 
Wind  and  yield  condition.  <j>.g(x)  was  integrated  trapezoidally 
as  before. 

The  eighteen  cases  examined  for  the  yield  and  Wind  con¬ 
ditions  specified  in  Table  I  conserved  activity  to  three  or 
four  significant  figures.  A  slight  reduction  in  recovered 
activity  was  noted  at  high  winds  for  each  yield  accompanying 
a  reduction  in  cumulative  <|>.g(x).  This  is  due  to  the  finite 
nature  of  the  downwind/upwind  integration  limits  signifying 
that  small  amounts  of  activity  still  remained  suspended  in 
the  cloud  at  the  completion  of  the  integration. 

To  complete  this  examination  of  conservation,  it  was 

also  necessary  to  verify  that  ,  /a<J>.g(x)dx  =  /ag(x)dx  in 

D  o 

Equation  (18)  for  the  yield  and  wind  conditions  used  in  this 
section.  To  do  so  required  setting  <J>  =  1.0  and  integrating 
Equation  (32)  as  before.  The  function  f  was  integrated 
analytically  from  -°°  to  +°°  in  the  manner  described  above. 

The  remaining  downwind  integration  was  again  accomplished 
numerically  using  the  trapezoidal  technique  described 
earlier  in  this  section.  This  time  however,  the  integration 
limits  were  from  ground  zero  (x  =  0)  to  the  .1  Roentgens/hr 
dose  rate  contour.  The  function  g(x)  was  separately  inte¬ 
grated  via  trapezoidal  integration.  Downwind  distance  (a) 
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TABLE  III 


Comparison  of  Distance  Traveled  Downwind  (a)  While 
Integrating  the  Fallout  Pattern  to  Recover  SNC  . 


Yield  for 

the  Present  Version 

of  f 

y 

Yield 

Wind 
(st.mi. ) 

(a)  with  g(x) 
(Nautical  Miles) 

(a)  with  4>.g(x) 
(Nautical  Miles) 

1KT 

0 

0.988 

0.957 

30 

120.6 

120.6 

60 

217.5 

217.5 

10KT 

0 

01.92 

01.87 

30 

485.6 

485.6 

60 

888.2 

888.2 

100KT 

0 

06.90 

06.62 

30 

1256. 

1256. 

60 

2330. 

2330. 

1MT 

0 

17.63 

17.23 

30 

2254. 

2254. 

60 

4237. 

4237. 

10MT 

0 

41.66 

40.79 

30 

3256. 

3256. 

60 

6180. 

6180. 

100MT 

0 

95.56 

93.75 

30 

4240. 

4240. 

60 

8104. 

8104. 

to  the  . 1  Roentegens/hr  dose  rate  contour  for  each  Wind 
and  yield  condition  was  compared  with  the  downwind  distance 
used  to  generate  the  data  in  Table  I.  The  step  sizes  used 
during  this  numerical  integration  were  identical  to  those 
used  earlier  for  each  yield  and  Wind  condition. 

In  the  eighteen  cases  examined  (specified  by  Table  1) 
the  recovered  activity  and  cumulative  g (x)  were  identical  to 
the  data  presented  in  Table  1  to  four  significant  figures. 
Also  as  Table  III  indicates,  the  downwind  limit  used  is 
identical  for  both  cases  to  four  significant  figures  for 
winds  greater  than  0.  For  a  zero  Wind  condition,  the  dispar¬ 
ity  fluctuates  between  3.1%  at  1  KT  in  "x"  to  1.9%  at  100  MT. 
This  is  not  significant  when  one  considers  the  integration 
technique  and  the  relationship  of  g(x)  and  <|>.g(x)  as  in 
Figure  (3).  Thus,  the  use  of  ./  <{).g(x)dx  or  /  g(x)dx  in 

D  0 

Equation  (18)  does  not  affect  model  conservation. 

Conclusion 

WSEG,  as  presented  by  Pugh  and  Galiano  in  their  origi¬ 
nal  work  is  mathematically  conservative.  It,  however,  from 
modifications  in  1962,  did  not  reflect  an  accurate  picture 
of  the  true  fallout  pattern  based  upon  later  data.  The  1962 
version  represents  this  data  more  adequately  but  as  demon¬ 
strated  in  this  section,  is  not  conservative  because  of  the 
addition  of  a2  to  f^.  The  effect  is  seen  primarily  at  yields 
less  than  1.0  MT  with  high  winds. 


V.  Computer  Implementation 

The  Fortran  subroutine  containing  the  WSEG  model  located 
in  Appendix  A  was  adapted  to  the  ASD  CYBER  74  computer  and 
coupled  with  a  main  program  designed  to  generate  output  iden¬ 
tical  to  the  sample  output  also  contained  in  Appendix  A  on 
pages  58-61.  All  work  was  subsequently  done  on  the  ASD 
CYBER  74  computer  and  specific  user  instructions  concerning 
program  operation  are  contained  in  Appendix  E.  Additionally 
this  program  is  designed  to  output  g(t),  <J>.g(t),  g(x),  cumu¬ 
lative  <f>.g(x),  the  g(t)  time  constant,  and  the  g(t)  exponent 
(n)  . 

The  results  are  contained  in  Appendix  B,  pages  72-75 
for  four  conditions  used  to  validate  the  computer  program. 

In  all  cases  the  yield  used  is  .01  KT  and  the  shear  is  .1 
knots/kilofeet .  The  effective  wind  varies  from  one  to  ten 
knots  and  the  fission  fraction  is  assumed  1.0.  As  seen,  the 
output  generated  for  this  thesis  is  nearly  identical  to  the 
output  in  Appendix  A  for  each  wind  condition.  Downwind  and 
crosswind  range  deviations  are  limited  to  .1  nautical  miles 
or  less  and  dose  rate  deviations  are  less  than  10  Roentgens/ 
hr. 

As  a  final  note,  the  output  created  during  this  inde¬ 
pendent  study  and  presented  in  Appendix  B  is  not  exactly 
identical  with  the  sample  output  of  Appendix  A  for  two  reasons 
One,  the  subprograms  programs  providing  the  cumulative  nor¬ 
mal  and  the  gamma  functions  were  locally  created  and  may  not 
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provide  the  same  accuracy  as  those  used  to  generate  the  sample 
output  in  Appendix  A;  two,  the  step  size  used  to  generate 
the  sample  output  was  unknown.  This  factor  plays  a  critical 
role  when  attempting  to  duplicate  earlier  work. 
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VI.  WSEG  Limitations 


The  purpose  of  this  section  is  to  discuss  several  limi¬ 
tations  of  the  WSEG  model  that  have  been  discovered  either  by 
researching  literature  or  through  experience  using  the  model. 
Some  of  the  following  limitations  have  been  mentioned  earlier 
in  this  report: 

1.  The  model  cannot  account  for  complex  wind  or 
shear  patterns.  This  restriction  leads  to  poor 
results.  Whether  results  are  high  or  low  de¬ 
pends  on  the  test  data  evaluated  and  the  accur¬ 
acy  of  the  wind  data. 

2.  The  model  is  unable  to  account  for  meteorologi¬ 
cal  conditions  such  as  rain,  snow,  etc. 

3.  Stabilized  cloud  parameters  appear  to  be  in¬ 
accurately  predicted  when  compared  with  other 
models  (Refs.  7:84;  and  5:17).  The  net  re¬ 
sult  is  a  reduced  fallout  pattern  area  as  both 
downwind  and  crosswind  displacement  are  affected. 
WSEG  also  underpredicts  and  cloud  center 
heights  at  low  and  high  yields.  This  inaccuracy 
is  partially  compensated  for  because  oq  is  over¬ 
predicted. 

4.  WSEG  produces  peak  concentrations,  for  the  con¬ 
dition  of  Wind  f  0,  that  occur  at  very  nearly 
the  same  downwind  location  regardless  of  wind 
velocity.  To  explain  this  error,  recall  that 
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the  downwind  distribution  function  is  <j>.g(x). 

The  parameter  "oti"  defined  by  Equation  (22)  with¬ 
in  <J>  varies  according  to  effective  wind  velocity. 
Higher  wind  velocities  lower  the  magnitude  of 
cti  thereby  driving  if  to  a  maximum  or  minimum 
more  rapidly  as  one  marches  downwind  or  upwind 
respectively  from  ground  zero.  Also  a  time 
limit  for  the  argument  of  <f>  was  established 
through  subroutine  logic  which  if  exceeded  set 


4>  =  1. 

Figure  (9)  depicts  this  discrepancy  clearly. 
It  represents  D^^/Yield  versus  distance  for 
several  wind  conditions  where  DH+^  is  the  hot¬ 


line  value.  Winds  are  in  units  of 


st .mi . 
hr 


The 


yield  condition  used  was  chosen  for  convenience 
during  preliminary  work.  Note  the  peaks,  when 
affected  by  effective  wind,  occur  at  very  nearly 
the  same  downwind  location.  In  fact,  the  peak 
for  the  60  st^1—  curve  even  appears  closer  to 
ground  zero  than  that  of  the  40  St^1:-  curve. 

This  error  indicates  that  either  the  ef¬ 
fective  wind  velocity  affects  the  fall  rates 
or  that  the  peaks  are  caused  by  large  particles 
whose  downwind  position  is  independent  of  wind. 
In  both  cases  the  inconsistencies  are  obvious. 

In  the  first  case  the  effective  wind  is  a  hori¬ 
zontal  contribution  to  the  translation  of 
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activity  and  has  no  effect  vertically.  The 
second  case  does  not  account  for  the  transition 
from  the  "0"  wind  condition  with  a  peak  at 
ground  zero  to  a  variety  of  wind  conditions  re¬ 
sulting  in  peaks  at  almost  the  same  downwind 
location. 


The  "0"  wind  curve  in  Figure  (9)  has  an  obvious 

asymmetry.  The  explanation  involves  the  torroi- 

dal  growth  expression  in  which  is  defined  as 
8|x  +  2a  | 


°o  (1  + 


-)  .  This  term  dominates  in  a 


"0"  shear  environment.  The  inclusion  of  2 o 

x 

within  the  absolute  value  sign  resulted  in  a 
reduction  in  the  magnitude  of  for  upwind 
calculations  when  compared  to  downwind  calcula¬ 
tion  at  the  same  |x|.  The  reason  is  that  ax 
is  always  positive.  This  effect  coupled  with 
the  three  hour  effectiveness  of  torroidal  growth 
produced  the  asymmetry.  A  correction  is  made 


simply  by  defining  the  torroidal  growth  expres- 
8  (|x|  +  2 a  ) 

sion  as  a  (1  +  - = - )  .  This  correction 

O  Li 

is  used  in  Figure  (10) .  Note  that  the  curves 
with  10,  20,  40,  60  —  winds  in  Figures 

(9)  and  (10)  are  identical.  This  discrepancy 


had  no  effect  in  downwind  fallout  pattern  com¬ 
putation  for  an  environment  with  wind. 


6.  The  use  of  the  Source  Normalization  Constant 
(SNC)  at  2  x  10£  Roe^tgens-mif,  adjusts  total 
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activity  deposited  as  local  fallout  at  80%  of 
fallout  entering  the  cloud  regardless  of  yield. 

7.  The  choice  of  Source  Normalization  Constant 
used  by  WSEG  is  three  to  four  times  greater 
than  used  by  later  models  (Ref.  7:85).  This, 
if  the  other  models  are  accurate,  would  cause 
an  overprediction. 

8.  The  model  generates  contour  patterns  that  are 
nearly  elliptical  in  shape.  As  seen,  the  com¬ 
puter  code  in  Appendix  B  solves  for  crosswind 
component  (y)  using  a  Gaussian  distribution. 
Complex  patterns  are  therefore  not  possible. 


48 


nit-  r 


VII .  Summary 


The  purpose  of  this  thesis  is  to  recreate  and  document, 
for  local  use,  the  most  popular  analytical  fallout  model, 
WSEG-10.  Additionally  several  sections  were  also  devoted  to 
analyzing  different  facets  of  WSEG-10.  This  study  will  pro¬ 
vide  a  basis  for  future  fallout  studies. 

The  first  section  discussed  the  WSEG-10  model  from  the 
original  document  by  Pugh  and  Galiano  (see  Ref.  1)  including 
later  revisions.  An  explanation  of  terms  was  provided  where 
possible. 

The  second  section  contained  an  analysis  of  the  cross¬ 
range  dispersion  term,  o^.  It  was  found  that  shear  effects 
predominate  at  late  times  after  burst  while  the  torroidal 
growth  term  is  dominant  soon  after  burst.  Graphs  of  several 
wind  and  shear  conditions  can  be  seen  in  Figures  (6) ,  (7)  , 

or  (8)  and  in  Appendix  C. 

Also  discussed  was  the  resemblance  of  the  torroidal  growth 
term  to  a  term  representing  diffusive  growth  based  on  Fick's 
Law.  Appendix  D  contains  a  development  of  Diffusivity  from 
Fick's  Law  for  comparison.  The  results  indicated  that  the 
process  defined  by  Pugh  and  Galiano  was  not  diffusive  growth 
for  the  following  reasons: 

i.  Arbitrary  deletion  of  units  in  the  torroidal  growth 
term.  Tcwas  assigned  a  dimensionless  value  of  8. 

This  produced  the  resulting  units  for  Diffusity. 
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2.  When  compared  with  diffusive  growth  modeled  by 
DELFIC,  the  WSEG  expression  for  torroidal  growth 
was  many  orders  of  magnitude  greater. 

3.  A  three  hour  limit  was  placed  on  the  effects  of 
this  term  thereby  restricting  its  contribution. 

The  third  section  discusses  the  property  of  conservation 
for  the  WSEG  model.  Results  demonstrated  that  activity  was 
conserved  regardless  of  the  upwind/downwind  normalized  dis¬ 
tribution  function  choosen  within  the  context  of  that  section. 


Only  with  the  original  crossrange  distribution,  f  ,  however, 
did  WSEG  conserve  activity.  The  1962  version,  which  is  also 
the  present  version  of  f^,  is  unnormalized  resulting  in  signi¬ 
ficant  losses  at  low  yields  and  high  winds.  It  was  also  found 
that  varying  shear  conditions  did  not  affect  conservation. 

The  fourth  section  describes  computer  implementation  of 
the  subroutine  Dose  obtained  from  Mr.  Ralph  Mason.  Subroutine 
Dose  contains  the  analytical  expressions  developed  in  the  WSEG 
model  consolidated  in  subroutine  form  for  easy  use.  Output 
nearly  identical  to  that  provided  in  Appendix  A  is  contained 
in  Appendix  B.  Appendix  B  also  contains  a  computer  listing 
of  the  AFIT  version  of  WSEG  along  with  a  definition  of  terms 
and  sample  results.  A  program  user's  guide  is  in  Appendix  E. 


The  last  section  discusses  a  series  of  weaknesses  or 


limitations  to  the  model  discovered  either  through  computer 
use  or  researching  the  literature  for  this  thesis.  Inconsis¬ 
tencies  covered  in  discussions  concerning  or  conservation 
were  not  included  in  this  section.  Several  of  these  weaknesses 
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included  a  D„,,  asymmetry  in  a  "0"  wind  condition  and  peak 

H“r  i 

D„  ,  at  nearly  the  same  downwind  location  for  winds  between 
H+ 1 


10  and  60 


st .mi . 
hr 
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Appendix  A 

Subroutine  Dose  and  Sample  Output 


This  appendix  contains  the  Fortran  subroutine  Dose  ob¬ 
tained  from  Mr.  Mason,  National  Military  Command  Support  Cen¬ 
ter.  Also  included  is  a  series  of  sample  results  for  the 
following  conditions : 

Yield  =  .01/. 03  MT 
Shear  =  .1  knots/kilofeet 
Fission  Fraction  =1.0 
Wind  -  1.0  to  10.0  KTS 

No  modification  has  been  done  to  the  subroutine.  References 
to  "dose"  in  the  sample  output  actually  refer  to  dose  rate 
in  Roentgens/hr. 
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oonnnooonno  n n o o n n  nooon  nnnnnnoooo  nnnnnnnoo  onoo 


SUBROUTINE  DOSE ( DB , DH , SIGYA2 , YY , XX , SHEARY , WIND , FFRAC , 
YIELD) 


THIS  SUBROUTINE  IS  THE  FEBRUARY  23,  1962  VERSION  WITH 
ALL  CHANGES  +  MODIFICATIONS  TO  RESEARCH  MEMORANDUM  10 
INCORPORATED.  USERS  WILL  BE  INFORMED  OF  ANY  LATER 
MODIFICATIONS. 

IN  NORMAL  FULL  CALLS, 

OUTPUT  PARAMETERS 

DB=THE  BIOLOGICAL  DOSE  IN  ROENTGENS  (INFINITE 
PLANE  DOSE) 

DH=THE  H+l  DOSE  RATE  IN  ROENTGENS  (INFINITE  PLANE 
DOSE) 

SIGYA2=THE  TERM  SIGMA  Y  SQUARED  IN  SQUARE  NAUTICAL 
MILES 

(SOMETIMES  USEFUL  IN  INTEGRATION  OF  DOSE  AREAS) 
INPUT  PARAMETERS 

YY=THE  CROSSWIND  DISTANCE  PERPENDICULAR  TO  THE 
WIND  DIRECTION  IN  NAUTICAL  MILES 
XX=THE  DISTANCE  ALONG  THE  X  AXIS  PARALLEL  TO  THE 
WIND  DIRECTION  IN  NAUTICAL  MILES.  (XX  IS  NEG¬ 
ATIVE  FOR  UPWIND  LOCATIONS) 

SHEAR=THE  CROSSWIND  COMPONENT  OF  SHEAR 
WIND=THE  EFFECTIVE  FALLOUT  WIND  IN  KNOTS 
FFRAC=THE  FISSION  FRACTION 
YIELD=THE  YIELD  IN  MEGATONS 

NOTE  THAT  CALCULATIONS  ARE  NOT  REPEATED  FOR  PARAMETERS 
THAT  HAVE  NOT  CHANGED.  THEREFORE,  THE  CALL  MAY  BE 
SHORTENED  TO  EXCLUDE  THOSE  PARAMETERS  AT  THE  END  OF 
THE  CALLING  SEQUENCE  THAT  REMAIN  THE  SAME. 

CALL  DOSE (DB,DH,SIGAY2, YY) 

THIS  SUBROUTINE  MAY  BE  USED  AS  A  FUNCTION  SUBROUTINE 
WITH  THE  VALUE  OF  THE  FUNCTION  EQUAL  TO  THE  BIOLOGICAL 
DOSE. 

ANSWER=DOSE ( DB , DH , S I GAY  2 , YY , XX , SHEAR , WIND , FFRAC , 
YIELD) 

IS  EFFECTIVELY  ANSWER=DB 

A  THIRD  USE  IS  TO  INPUT  XX  AND  THE  DOSE  AND  RECEIVE 
AS  OUTPUT  THE  CORRESPONDING  YY  IN  NAUTICAL  MILES. 
(USEFUL  IN  COMPUTATION  OF  FALLOUT  CONTOURS) 

CALL  DOSE (YDH , YDB , -DOSE , YY , XX , SHEAR , WIND , FFRAC , 
YIELD) 

YDH=THE  YY  DISTANCE  IN  NAUTICAL  MILES  FOR  AN  H+l 
INPUT  DOSE 

YDB=THE  YY  DISTANCE  IN  NAUTICAL  MILES  FOR  A  BIO¬ 
LOGICAL  DOSE 

-DOSE=MINUS  THE  VALUE  OF  THE  DOSE 

ALL  OTHER  PARAMETERS  ARE  THE  SAME  AS  ABOVE 
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IF (YIELD-OLDYLD)  1,2,1 

C  YIELD  DEPENDENT  CALCULATIONS 

1  OLDYLD=YIELD 
YMT=LOGF (YIELD) 

T3=2000000. *YIELD 

SIGO=. 7+YMT/3.-3. 25/ (4 . * (YMT+5.4) **2) 

SIGO=EXPF (SIGO) 

SIG02=SIGO*SIGO 
Tl=YMT+2 . 42 

H=44.+6.1*YMT=. 205*Tl*ABSF (Tl) 

SIGH-. 18*H 
SIGN2=SIGH*SIGH 
T2=H/ 60 

T= (12 . *T2-2 . 5*T2*T2) * (1 . - . 5*EXPF (- (H/25. ) **2) ) * 
1.0573203 
GO  TO  3 

2  IF (WIND-OLDWIND)  3,5,3 

C  WIND  DEPENDENT  CALCULATIONS 

3  OLDWIND=WIND 

4  ZL0=WIND*T*1. 151515 
ZL02=ZL0*ZL0 

SIGX2=SIG02* ( ZL02+8 . *SIG02) / ( ZL02+2 . *SIG02) 
SIGX=SQRTE(SIGX2) 

ZL2=ZL02+2.*SIGX2 
ZL=SQRTE ( ZL2 ) 

T14=ZL02+.5*SIGX2 
40  ZN= (ZL02+SIGX2J/T14 

IF (ZN-1. 002)  102,102,103 

102  ZN=1. 

T20=l. 

GO  TO  42 

103  T20=GAMMA(1.+1./ZN) 

42  T4=T3/ (ZL*T20*2 . 5063) 

PALPH=. 001*H*WIND*1 . 151515/SIGO 
ALPH1=1 . / ( 1 . +PALPH) 

T5=ZLO/(ZL*ALPHl*SIGX) 

T6=2. *SIGX2*T*T*SIGN2/ZL2 
T15=ZL02/ZL2 
T7=T15*T*T*SIGH2 
GO  TO  6 

5  IF (SHEARY-OLDSHR)  6,8,6 

C  SHEAR  DEPENDENT  CALCULATIONS 

6  OLDSHR=SHEARY 

7  T21=SHEARY*SHEARY*1. 325975 
T8=T6*T21 
T9=T7*T21/ZL2 

GO  TO  9 

8  IF (XX-OLDX)  9,116,9 

C  X  DEPENDENT  CALCULATIONS 

9  OLDX=XX 
X=XX+6080 ./5280 . 

10  Tl0=X+2. *SIGX 

Tll=l.+ (8. *ABSF (T10) )/ZL 
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IF (Tll-4 . )  12,12,11 

11  Tll=4 

12  T22=Tll*SIG02 
T30=T5*X 

IF (T30-6 . )  35,36,36 

36  T30=l . 

GO  TO  37 

35  T30=CUMNOR(T30) 

37  T12=T9*T10*T10 

43  IF ( X)  13,14,13 

14  T13=1 . 

TO  TO  15 

13  IF ( ZN-1 . )  113,114,113 

114  T13=EXPF (- (ABSF (X) /ZL) ) 

TO  TO  15 

113  T13-EXPF (- (ABSF (X) /ZL) **ZN) 

15  SIGY2=T22+T8+T12 
SIGY=SQRTE (SIGY2) 

TARR=SQRTE ( .  25+ (T15*T10*T10*T*T*2 . *SIGX2)/T14) 
BETA=LOGE (TARR/31 .6) 

ZLD=- . 287- • 52*BETA- . 04475*BETA*BETA 
BIO=EXPE (ZLD) 

IF (WIND)  27,27,53 

53  T23=(2.*X)/(WIND*1. 151515) 

IF (T23-10 . )  28,28,27 

27  ALPH22=1. 

GO  TO  29 

28  T24=CUMNOR(T23) 

ALPH22=1 . / ( 1 . +PALPH* ( 1 . -T24 ) ) 

29  ALPH2=ALPH22*ALPH22 

IF (SIGYA2 )  91,90,90 

TO  CALCULATE  V,  GIVEN  X  AND  A  DOSE 

90  SIGYA2=SIGY2*ALPH2 
GO  TO  17 

91  DHX0=T30*T13*T4*FFRAC/SIGY 
DBXO=DHXO*BIO 
DOSEL=ABSF ( S IGYA2 ) 

IF (DOSEL)  95,117,95 

95  IF (DHXO/DOSEL-1 . )  94,94,92 

94  DH=0. 

GO  TO  93 
117  DH=0. 

DB=0. 

RETURN 

92  DH=ALPH22*SIGY*SQRTE(2. *LOGE (DHXO/DOSEL) ) *5280. /6080 

93  DB=ALPH22*SIGY*SQRTE (2 . *LOGE (DBXO/ DOSEL) ) *4380 ./6080 
RETURN 

116  IF (SIGYA2)  91,16,16 

16  IF ( YY-OLDY)  17,19,17 

C  Y  DEPENDENT  CALCULATIONS 

17  OLDY=YY 
Y=YY*6080./5280. 

18  T16=EXPE(-.5*Y*Y/(ALPH2*SIGY2) )/SIGY 
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20,22,20 


GO  TO  22 

19  IF (FFRAC-OLDFRAC) 

20  OLDFRAC=FFRAC 

21  DH=DDH*FFRAC 
DB= DDB  *FF  RAC 
RETURN 

22  DDH=T30*T16*T13*T4 
DDB=DDH*BIO 

GO  TO  20 

END 

END 
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Appendix  B 

AEIT/WSEG  Fortran  Computer  Program 

This  appendix  contains  the  fully  documented  AFIT/WSEG 
Fortran  computer  code  utilizing  the  subroutine  Dose.  It  also 
contains  a  definition  of  terms  for  Dose  and  sample  output 
for  four  different  conditions.  Disregard  the  computer  gen¬ 
erated  sequencing  at  the  left  margin  as  it  is  not  essential 
to  the  operation  of  the  computer  code. 


1H‘ 
ltd* 
120- 
IN- 
144* 
ISO- 
tfP-.' 
179-C 
IM-C 
J90-C 
2M-C 
2l#-C 
22#-C 
4N-C 
24«-C 
254-C 
2W-C 
27S-C 
2M-C 
296-C 
3M*C 
31#-C 
320 -C 
330 -C 
340-C 
3S#-C 
360 -C 
37# -C 
3K0-C 
390-C 


PROGRAM  WSEG<  INPUT*/##. OUTPUT. TAPIS-OUTPUT) 

DP1ENSI0N  D4*l  2#  J.UPNAXI  (2# ) .  DUPHAX1  ( 2#  ).  YYNAXI  <2# ) 

DIWMSION  R2NAXU1 (2# l.CUNCXI <2#  > 

COHflON  GT.CX.CT1 

common  oldvld.oldwind.oldshr.oldx.oldy, oldfrac 

CONHOH  T.2I4 

TMI*  PfNV.OAP  IS  THE  AIR  FORCE  INSTITUTE  OF  TECHNOLOGY. S  UERSIOH 
OF  WStG-1#  CREATES  »V  PUON  IN  IBS 9  AND  CONDENSED  INTO  PRESENT  FORM 
--TFC  SUIROUT INE  DOSC--OSTAINED  f  RON  NR.  RALPH  NASON  (NATIONAL 
NILITARV  CONNAND  SUPPORT  CENTER).  INPUT  AND  OUTPUT  PAR  (PETERS  TO 
DOSE  ARE  EXPLAINED  UITHIN  DOSE  AND  NOT  REPEATED  HERE.  ADDITIONAL 
C04EENTS  HAKE  SEEN  ADDED  WITHIN  DOSE  TO  TO  FURTHER  AID  THE  USER  IN 
IDENTIFYING  THE  VARIABLES  AND/OR  THE  PROCESS  1NU0LUED. 

TtC  NAIN  PROGRAM  IS  DIUIDCD  INTO  THREE  SECTIONS) 

1.  THE  FIRST  DO  LOOP  SEARCHES  D0UNU1ND  OF  GROUND  ZERO  FOR 

NAY I NUN  UNIT  TINE  REFERENCE  DOSE  RATE  ECT.  IT  ALSO  CALCULATES 
CUNULAT1UE  G(X>  . 

2.  THE  SECOND  DO  LOOP  SEARCHES  UPWIND  OF  GROUND  ZERO 

3.  THE  THIRD  DO  LOOP  SEARCHES  FOR  CROSSRANQE  DATA.  IT  IS 

ONLY  CONCERNED  WITH  DOwNUIND  DUE  TO  THE  NATURE  OF  THIS  NOOEL. 


INPUTS  (StSlttt>t««(Sf*KtS(*«fS*tSUttSt*tS»S(**»Stt>SSS(St(»SaS«SS 

FFRAC- -REAL  NUNDER  SPECIFYING  FISSION  FRACTION  FOR  DURST 

I  YIELD- -INTEGER  PARAMETER  SPECIFYING  THE  NUNKR  OF  YIELDS  TO  K 
EUALUATED. 


41#-C  I  SHEAR- INTEGER  SPECIFYING  THE  NUNKR  OF  SHEAR  CONDITIONS. 

42*>-C 

430-C  IWIND--  INTEGER  SPECIVINO  THE  NURSER  OF  WIND  CONDITIONS 

44#*C 

4SR-C  IGT-- INTEGER  SPECIFYING  OUTPUT  1NCLUD1N0  G(  T )  AND  TINE.  IF 
46# -C  DESIRED  ENTER  I,  IF  NOT  ENTER  • 

47R-C 

4B0-C  I GX-  INTEGER  REQUESTING  OUTPUT  OF  G(X)  AND  DOWNWIND  DISTANCE.  IF 

4N-C  DESIRED  ENTER  i.  IF  NOT  ENTER  t 

S##*C 

St«-C  ICUNGX—  INTEGER  REQUESTING  CUNULATIUC  C(X)  FOR  EACH  I (R-UT  DOSE 

52# -C  RATED  CONDITION.  IF  DESIRED  ENTER  I.  IF  NOT  ENTER  S 

53#-C 

S4#-C  XLEN — REAL  NUNSER  SPECIFYING  THE  DOUHWIM  AND  UPWIND  MARCH  INC 

SS#*C  INTERVAL.  THE  UNITS  ARC  NAUTICAL  NILES. 

ss#-c 

S7#-C  IHT- -INTEGER  SPECIFYING  WHICH  ITERATION  THE  WRITE  STATEMENTS  FOR 

S##*C  G  ( X )  AND  «(T)  ACT  (TON.  I.E.  IF  INT-l#  THEN  EVERY  TENTH  VALUE  OF 

S##*C  G(X)/G(T)  AND  DISTANCE /TINE  WILL  K  PRINTED. 

G##*C 

S1#*C  YIELD- -NEAL  NUNRCR  SPECIFYING  THE  YIELD  OF  THE  WEAPON  IN  NEQATONS. 
S2#*C 

S3#-C  WIND — REAL  NUNKR  SPECIFYING  THE  EFFECTIVE  WIND  IN  KNOTS. 


tH-C  SHE  ARY-  -REAL  NUNSER  SPECIFYING  TK  CROSSWINS  SKAR  COFTOKWT  IN 

(M-C  XNOTS/TCILOFOOY . 

S7S-C 

CM-C  DHI — REAL  NUNKR  SPECIFYING  THE  UNIT  TIK  REFERENCE  DOSE  RATE 

GOS-C  THE  COMPUTER  WILL  USE  AS  IT  GENERATES  T«  OUTPUT  PARAfCTCM. 

7s#»e 
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710-C 

720  *C 

730-c  >:MWitiiiuuiuiiut'.<<utiittitmuiitiiiiiimit>iiii«iiimii 

7  40-C 

??0*C  INITIhL  cohditions—vield.uind, shear. fission  fraction, step  size 

7sc*.: 

— *'•(-  V.  i--0€POf  ITION  OF  FALLOUT  PER  LINEAR  NILE.  THE  UNITS  ARE  PER 

730 -C  NAUTICAL  NILE.  INCLUDED  IS  CORRESPONDING  RANGE  FROM  GROUND 
7A>*0  ZERO. 

5W*C 

310-C  G<T>  — DEPOSITION  OF  FALLOUT  EUERYUCRC  PER  TINE.  THE  UNITS  ARE 
320*C  PER  HOUR.  ACCONPANV1NG  ClT )  IS  ITS  TINE  COORDINATE  IN  HOURS. 

S30-C 

340-C  DUDNAX- -DISTANCE  IN  NAUTICAL  NILES  FRON  GROUND  ZERO  DOIJNUIND  TO 

SS0-C  UTRD  RATE  SPECIFIED  IV  DNI . 

860 -C 

870 -C  UPNAX- -UPWIND  DISTANCE  TO  UTRD  RATE  SCECIF1ED  IV  DHI  FRON  GROUND 

380-C  ZERO.  THE  UNITS  ARE  NAUTICAL  NILES.  NOTE  XX  IMV  IE  <-)  OR  (♦> 

S90-C  DEPENDING  ON  THE  NAGNITUDES  OF  THE  EFFECTIVE  UIND  AND  YIELD. 

900 -C 

910-C  DINAX — NAX1NUN  HOTLINE  UTRD  RATE  CONTAINED  UITMIN  THE  TOTAL 

920«c  fallout  pattern  specified  iy  the  nininun  dhi.  units  we  r  hour. 
930*C 

9-tO-C  RHAXD— DISTANCE  FRON  GROUND  ZERO  IN  NAUTICAL  NILES  TO  DDNAX. 

950  *C 

960-C  DGZ— UTRD  RATE  AT  GROUND  ZERO.  THE  UNITS  ARE  R/HOUR. 

970 -C 

980-C  VVHAX— NAX1NUN  CROSSRANGE  UIDTH  OF  ISO-DOSE  RATE  CONTOUR  SPECIFIED 
990-C  »v  DNI.  THE  UNITS  ARE  NAUTICAL  NILES. 

1 000-C 

1010-C  R2NAXU — DOUHUIND  OR  UPWIND  DISTANCE  TO  YYTWX.  UNITS  ARE  NAUTICAL 

I02O-C  NILES. 

1030-C 

1040-C  CUNCX— CUHULATIUE  G(X>  SOUNDED  DY  TX  UPWIND  AND  DOWNWIND  RANGE 

iese-c  data  specified  iy  dhi.  cuncx  calculated  dv  trapezoidal  intec- 
1060-C  RATION  and  is  dinens ionless. 

1070-C 

1080-C  -N-— EXPONENT  OF  C<X>  OR  CITI 

1090-C 

1100-C  ’T’— TINE  CONSTANT  FOR  YIELD 

U10-C 

1120-C 

1130-C 

1140*C 

1150*  RCAOt.FFRAC 

USA*  SICYA2*0.0 

1170-  DO  10  J-1.20 

1181-  DHRU>*0.0  SUPNAXKJ  1*0.0  •  YYNAXK  J 

1190*  R2NAXUK  J)-0.0  •  DUDNAXKJ  1*0.0  0  CUNOXI  fj  )*0.0 

1200*10  CONTINUE 

1210*  READ* .  I  YIELD.  I  SHEAR,  IUIND.  IGT,  ICX,  ICtMX 

1220*  READ*  ,XL£H,  1NT 

1230-  DO  14  JK*1 ,  I  YIELD 

1240*  READ*. YIELD 

12S0*  DO  14  JKIC  ■  1 ,  I  SHEAR 

1280*  READ*,  SMEARY 

1270*  DO  14  ICJK*I, IUIND 

1280*  READ*. UIND 

1290*  WRITE  2S 
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1300-25 

1310- 

1320-30 

133.'* 

1340-35 

1330- 
1 ?Ov *50 

1370- 

13'0-SS 

1350- 

1400-60 

141C- 

1420*64 

1430- 

1440-66 

1450- 

1460-66 

1470- 

1480-67 

1490- 

1500- 

1510- 

1620-68 

1S30- 

1640-69 

1550- 

1560- 

1570- 

1530- 

1590- 

1600- 

1610- 

1620- 

1630- 

1640- 

1650- 

1660* 

1670- 

1680- 

1690- 

1700- 

1710- 

1720- 

1730- 

1740- 

1750- 

1760- 

1770- 

1780- 

1790- 

1800- 

1810- 

1820- 

1830- 

1840* 

I860- 

I860- 

1870- 

1880- 

1890- 


FORHAT» 1X.*CALCULATED  M4l  HOUR  DOSE  RATE  CONTOURS*.// > 

UR1TE  30 

FORHAT! IX, I  INITIAL  CONDITIONS)*./ » 

UUTt  35. YIELD 

r0RHATiS<.mELD<HEGAT0NS>*.18X,»-  t,F6.2,/> 

UPITE  50.FFRAC 

F.-’RIViT(5X,»FISSI0N  FRACTION*. 9X,t-  ( ,F6.2.  /  > 

UR1TE  55.U1N0 

FORMAT (5X,*UIND!KTS1*.1€X,*-  X.F6.2./) 

URITE  60,  SHEARV 

F OPMAT !5X, (SHEAR! KTS  PER  K1LOFOOT  >t.2X,(-  *.F6.2,/I 
UR1TE  64.XLEN 

FORMAT !SX,*ST£P  SIZE1NAUTICAL  NILE)*,*  •  *,FS.3,/| 

URITE  65 

FORMAT! IX. XRESULTS)  ALL  DISTANCES  IN  NAUTICAL  NILES*. //> 

IF' IGT.EO. 1)  URITE  66 

FORMAT!  10X.*G(T)-.1SX,  *01  T)*PHI*.14X, ‘TINE*  1 
IF i IGT.EO. 1 >  URITE  67 

FORMAT (7X.XPER  HOUR*.  14X. *PCR  HOUR* .  14X, (HOURS*,//  ) 
IFlUIHO.EO.O.e.ANO.IGT.EO.l  )  PRINT*. *G(T)  AND  TINE  ARC  FUNCTIONS 
'OT  U1NP  AND  DISTANCE  AND  ME  EITHER  UNDEFINED  OR  8.8* 

IFllCX.EO.ll  URITE  68 

FORMAT! 7X.*C!X  )*.  1SX,  (HOTLINE*  ) 

IFlIGX.EO.ll  URITE  69 

FORMAT I5X.»PER  NAUT.N1  .  *  ,6X,  (FROM  CRD. ZERO*. // ) 

D8MAX-0.0 

DO  4  JJ-1.8 

CUNGX-8.0 

UPMAX-0.8 

XX-8.0 

VY-0.0 

OLDFRAC  -0.8 

0LDY-1.E9 

0LDX-1.E9 

0LDSHR-1.E9 

0LDUIND-I.E9 

OLDYLD-8 .8 

DM0LD-0.8 

DH-I.E9 

READ*. DM1 

DO  1  1-1.6000 

KL-1I/INTUIHT 

CALL  DOSE ( DO .  DH . S I G V A2 , W , XX . SHEARV , UIND.FTRAC .YIELD  1 

IF!YY.E0.8. .AND. XX. EG. 8. 1  DCZ-DH 

1F<DN.LT.DNI.AND.I.E0.1  1  DUDMAX-8.8 

IF  I  DM. LT.DN1 .AND. I .EG. 1  1  YVNAX-8.8 

lFIDH.LT.DHI.AND.I.EO.l  I  R2NAKV-8.G 

1FIDH.LT.DNI. AND. I. EG. 1  1  CUNGX-8.8 

IF IDH.GE . DINAX 1  DIMAX-DH 

IF!  DH.GE  .  D8NAX  1  RMAXD-XX 

IF(  DH.CE.DHI ,  AND.  DH .  LE  .  DHOLD  I  DWDNAX-XX 

IF!DH.LE.DN1.AND.DN.CC.DH0LD.AND.D8NAX.GC.DHI  I  UPNAX-XX 

IFIDH.LT. DNI  .AND.DN.LT.  DHOLD  I  GO  TO  6 

IFII.GT.1I  AU-ICXOLD-GMI/3. 

IFiI.GT.il  CUNCX-AUSXLEN4CUN6X 
CXOLD-CX 

IF IUIND.EO.O.O I  TIME -8. 8 
IF iUIND.GT. 8.8  >  TINE-XX/U1ND 
IFIJJ.CT.il  GO  TO  9 

IFiWIND. NE. 9. 9. AMD. IGT.EO. I. AND.KL.CO.il  URITE  78.GT.GTI.TINE 
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,  1 


1900-70  FORPIAT(3X,F12.8,10X,F12.8.12X,F6.2) 

1910*  IF<Kl.EQ.I.AND.ICX.E0.1 )  URITE  71. CX.XX 

1920-71  FORnAT<3X.F12.8,10X,F6.2> 

1930-9  IF  I  DFflAX .  LT  .  DHI  )  GO  TO  5 

1940-  PHOLD-DH 

19FC-  x;:-xx*xlen 

1960-1  CONTINUE 

1970-  IF<DH.GT.DHl.AND.I.E0.bee0)  URITE  72, DHI 

19S0-72  FORNATd0X,F6.0, 10X,tCALCULATI0NS  INCOMPLETE— PLEASE  INCREASE 

1990.  'STEP  SI2E*./) 

2000-  1F'DH.GT.DHI.AND.I.EQ.6000)  GO  TO  15 

2010-5  XX-0.0 

2020-  W-0.0 

2030-  IFdGX.EO.l.AND.JJ.EO.l)  PRINT*,-  • 

2040-  IFdGX.EO.l.AND.JJ.EO.l)  PRINT*.*IN  THE  UPUIND  DIRECTION!- 

2050-  IFdGX.EQ.l. AND.JJ.EO.il  PRINT*,-  ■ 

2060-  OLDYLD-0.0 

2070-  0LDUIND-1.E9 

2080-  0LDSHR-1.E9 

2090-  0LDX-1.E9 

2ie0-  0LDY-1.E9 

2110-  OLDFRAC-0.0 

2120-  DH-DBNAX 

2130-  DO  2  1-1,6000 

2140-  KL-d/IHTItlNT 

2150-  IF(DH.LT.DHI.AND.I.GE.21  GO  TO  7 

2160-  CALL  D0SE(D8,DH.SIGYA2.YY,XX.SHEARV.UIND.FFRAC, YIELD) 

2170-  IFtDH.CE.DHl)  UPHAX-XX 

2180-  IF(  DH.  LT .  DHI  .AND.  I  .EQ.  1 .  AND.UPflAX.EO.0.0  1  UPMAX-0.0 

2190-  IF(DH.GE.DBflAX)  DBNAX-DH 

2280-  IF(DH.CE.DBflAX)  RTIAXD-XX 

2210-  IFd.GT.ll  AU-  (GXOLD+GX 1/2, 

2220-  IFd.GT.ll  CunCX-AM*XLEN+COnGX 

2230-  CXOLD-GX 

2240-  IFCJJ.CT.l)  GO  TO  C 

22S0-  IFIKI.EO. I. AND.IGX.EO.il  URITE  71.GX.XX 

2260-6  XX-XX-XLEN 

2270-2  CONTINUE 

2280-7  XX-0.0 

2290-  VY-0.0 

2300-  OLDFRAC-0.0 

2310-  YYOLD-0.0 

2320-  0LDV-1.E9 

2330-  0LDX-1.E9 

2340-  0LDSHR-1.E9 

2350-  0LDUIND-1.E9 

2360-  OLDYLD-0.0 

2370-  DHII--DHI 

2380-  DO  3  1-1,6000 

2390-  IF1XX.GT  .DUDflAXl  GO  TO  11 

2400-  CALL  DOSEIYDB.YDH.DHII.YY, XX, SHEARY.UIND.FTRAC, YIELD) 

2410-  IF(  YDH.GE.  YYOLD  )  YYHAX-YDH 

2420-  IFtVDH.GE. YYOLD)  R2HAXU-XX 

2430-  YYOLD -VDH 

2440-  XX-XX4XLEN 

24S0-3  CONTINUE 

2460-11  R2DAXUI  ( JJ  1-R2NAXU  0CUT1GXI  ( JJ  1-CUflCX 

2470-  DUDHAXI  ( JJ  1-DUDNAX  S  VYT1AXI  ( JJ  l-YVHAX 

2480-  DHR<  JJ  l-DHI  >UP)MXI(JJ  1-UPNAX 

2490-4  CONTINUE 
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2530 
2510 
£=20 
2=20 
2"  4-: 
£550 
2=60 

£533 

2S53 

2600 

2610 

2620 

2620 

2640 

2650 

2660 

2670 

2630 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2730 

2790 

2200 

2810 

2820 

2830 

2840 

2250 

2860 

2870 

2380 

2290 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 


3060 

3070 

3080 

3090 


76 

13 


SO 


81 


85 

90 

95 

100 

185 

14 

15 


IF dCUMCX.EQ .  1 >  WRITE  73 

F0PMAT<  ",10X,*D0SE  RATE*, 1 OX , *M AX I  MUM*.  10X,*MAX1MUNX,  10X, XMAXINUN 
1  (.10  :,XPAHGE  TO*,10X,*CUnULATlUE*) 

IF-  ICUMCX.EO.l *  CO  TO  12 
UFITE  74 

fopmmTc  ,iox,*pose  rate*. iox.xmaxinumx.iox, maximum*. lex.xNAxinun 

'*.10N',*PANGE  TO*  1 
IF- lCIWJX.EO.ll  WRITE  75 

FOPHAT f  SX.  *R0£MTGENS/HRX,9X,XUPWINDX,  11X,  SD0WNUINDX.9X, XCROSSUINDX 
*.3X,XMA<  WIDTH*,  12X,*G(X>*,/> 

IF* ICUMCX.EO.l)  GO  TO  13 
UPITE  76 

FORMATS 8X,  *R0CNTGEM5/MR*,9X,*UPWIND*,  11X,SD0UNU2ND*,9X,  »CROSSUINDS 
•,8X,*HAX  WIDTH*, /I 
DO  8  01,8 

IFUCUMGX.NE.l)  WRITE  80, DHR(K l.UPNAXI  IK  ),  DUDMAXHK ),  YYMAXI  IK 
!  >,R2MAXUHK) 

FOPttATl  10X,F6. 0, 13X,F6.2,11X,FE,2, 1 1X,F6.2, 1HX.F6.2,/  ) 
IFdCUnGX.EO.il  UR1TE  81 .  DHR(K  )  .UPflAXl  (K  ) ,  DUDMAXI  (K  ),  YYTIAXI  (K  1 

•  ,R2HAXUI(IC),CUMCXI<K) 

FORMAT  (1OX.F6.0, 13X.F6.2, 11X.F6.2, 11X.F6.2, 12X.F6.2, 12X.F8.6, 

•  / ) 

CONTINUE 
PRINT*.*  * 

URITE  85.DBMAX 

format<sx, maximum  dose  rate*.hx.*-*,2x,f6.o,/) 

URITE  90.RMAXD 

FORMAT  ( 5X ,  tRANGE  TO  MAXIMUM  DOSE  RATE  •  *,FS.l,/> 

URITE  95.DGZ 

FORMAT < 5X , XDOSE  RATE  AT  CROUND  EERO  ■  X.F6.D,/) 

URITE  100, 2N 

F0PNAT<5X, XPARAMETER  *N**,lSX,*-*.2X,F6.4,/> 

URITE  105, T 

FORMAT (SX.XTinE  CONSTANT*. 15X,*-*, 2X,F7. 4, /// ) 

CONTINUE 

STOP  ‘END  OF  PROCRAM* 

END 

SUBROUTINE  D0SE(DB,DH,SIGYA2,  YY,  XX,  SHE  ARY,  UIND.FFRAC,  YIELD) 

COMMON  CT.GX.GT1 

COMflON  OLD YLD . OLDU I ND , OLDSHR , OLDX.OLDV, OLDFRAC 
COMMON  T.ZN 

THIS  SUBROUTINE  IS  THE  FEBRUARY  23,1962  VERSION  UITH  ALL  CHANGES/’ 
MODIFICATIONS  TO  RESEARCH  MEMORANDUM  19  INCORPORATED.  USERS  UILl 
BE  INFORMED  OF  ANY  LATER  MODIFICATIONS. 

IN  NORMAL  FULL  CALLS. 

OUTPUT  PARAMETERS  - 

DB-THE  BIOLOGICAL  DOSE  IN  ROENTGENS  FOR  AN  INFINITE  PLANE  DOSE. 

DH'THE  H+l  HOUR  DOSE  RATE  IN  ROE  NT  GENS /HOUR  FOR  AN  INFINITE  PLA 
DOSE. 

SIGYA2-THE  TERM  SIGMA  V  SQUARED  IN  SQUARE  NAUTICAL  RILES. 
(SOMETIMES  USEFUL  IN  INTEGRATION  OF  DOSE  AREAS) 


INPUT  PARAMETERS 
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31^0-C 

3110-C 

3120'C 

3130-C 

3H0*C 

3150-C 

3160-C 

3170-C 

31S0-C 

3190-C 

3200-C 

3210-C 

3220 *C 

3230-C 

3240-C 

32S0-C 

3260*C 

3270 -C 

3230-C 

3290-C 

3300-C 

3310-C 

3320-C 

3330-C 

3340-C 

33S0-C 

33G0-C 

3370-C 

3380-C 

3390-C 

3400-C 

3410-C 

3420-C 

3430-C 

3440-C 

3450-C 

3460-C 

3470-C 

3480-C 

3490-C 

3S00-C 

3510-C 

3S2C-C 

3S30-C 

3S40-C 

3550- 

3560*1 

3570* 

3580* 

3590* 

3600-C 

3610- 

3620* 

3630* 

3640-C 

3650* 

3660*C 

3670* 

3680* 

369#* 


YY-THE  CROSSUIND  DISTANCE  PERPENDICULAR  TO  THE  UIHD  DIRECTION 
IN  NAUTICAL  NILES. 

XX-PISTANCE  ALONG  THE  HOTLINE  PARALLEL  TO  THE  UIND  DIRECTION 
IN  NAUTICAL  NILES.  (XX  IS  NEGATIVE  FOR  UPWIND  LOCATIONS) 

SHE ARY* THE  CROSSUIND  COMPONENT  OF  SHEAR. 

UIMD'THE  EFFECTIVE  FALLOUT  UIND  IN  KNOTS. 

FFRAC-THE  FISSION  FRACTION. 

VIELfi-THE  YIELD  IN  MEGATONS. 


NOTE  THAT  CALCULATIONS  ARE  NOT  REPEATED  FOR  PARAMETERS  THAT  HAVE 
NOT  CHANCED.  THEREFORE.  THE  CALL  MAY  BE  SHORTENED  TO  EXCLUDE 
THOSE  PARAMETERS  AT  THE  END  OF  THE  CALLING  SEOUENCE  THAT  REMAIN 
THE  SANE. 

FOR  EXAMPLE*  CALL  D0SEtDB.DH.SICVA2.YY) 


A  SECOND  USE  IS  TO  INPUT  XX  AND  UTRD  RATE  AND  RECEIVE  AS  OUTPUT 
CORRESPONDING  YY  IN  NAUTICAL  MILES.  (USEFUL  IN  COMPUTATION  OF 
FALLOUT  CONTOURS). 

FOR  EXAMPLE  I 

CALL  DOSE<  YDH.VDB, -DOSE,  YV,  XX,  SHE  ARY.  UIND  FFRAC, YIELD) 
YDH-THE  YY  DISTANCE  IN  NAUTICAL  MILES  FOR  AN  H«1  INPUT  DOSE. 
YDB-THE  YY  DISTANCE  IN  NAUTICAL  MILES  FOR  A  BIOLOGICAL  DOSE. 
-DOSE-MINUS  THE  VALUE  OF  THE  DOSE. 

ALL  OTHER  PARAMETERS  ARE  THE  SAME  AS  ABOVE. 


NOTE  NAUTICAL  MILES  ARE  CONVERTED  TO  STATUTE  MILES  IN  DOSE 
THE  CONVERSION  FACTOR  -  1.151S1S  STATUTE  MILES  PER  NAUTICAL  MILE 


yield  dependent  eouations  ummiutmmumtm 
IF(YIELD-OU)YLD)  1.2.1 
OLDVLD* YIELD 
YNT  *ALOG( YIELD ) 

T3-2000000. XYIELD 

SI  GO*  .  7+Yrnv3 . -3 . 2S/ ( 4 .  ♦  C  YHT*5 . 4  >**2  ) 

INITIAL  STABILIZED  CLOUD  RADIUS 
SIGO-EXP(SICO) 

SIC02-SIC0TSIG0 

Tl-YNT+2.42 

CLOUD  CENTER  HEIGHT  IN  KILOFEET 
H-44 .  *6 . 1 IYNT- .  20STT 1 TABS ( TI  ) 

VERTICLE  STANDARD  DEVIATION 
SIGH*. 1804 
SIGH2*  SIGHTS  IGH 
T2-H/60. 


j 


68 


3700-C  TIME  CONSTANT 

3710-  T-(ia.»T2-2.S*T2*T2>*(l.-.S*EXP(-(H/25.  UJ2))*l.eS73203 

3729*  GO  TO  3 

3730-2  lFlWIND-OLDUIND )  3,5,3 

3?40*C  WIND  DEPENDENT  CALCULATIONS  ****»»XX*X***XX*X*t*****« 

37S0-3  OLDUIND-UIND 

3760-4  ZL0-UINDXTX1 .151515 

3770*  ZLC2-ZL0XZL0 

37S0*  SIGX2-S IG02X( ZL02-8 . XSXG02  >/t  ZL02+2.XSIG02 ) 

37S0-C  STANDARD  DEUIATION  FOR  X  DIRECTION 
3800-  SIGX-SIGX2XX.5 

3810*  ZL2-ZL02+2.XSIGX2 

3320-C  PARAMETER  “L* 

3830*  ZL*ZL2Xt.S 

3840*  T14-ZL02* .5XSIGX2 

3350"C  THE  PARAMETER  *N* 

3860-40  ZN-(ZL02*SICX2J/T14 

3870*  IF (ZN-1 , 002 )  102,102,103 

3880-102  ZN-1. 

3890-  T20-1. 

39O0-  GO  TO  42 

3910-C  GAHMA<1.+1./ZN)  FOLLOUSl 

3920-103  T20-GAflflAtl./ZN) 

3930-42  T4-T3^(ZLXT20X2. 50663) 

3340*  PALPH- .O01XHXUINDX1 . 1S1S15/SIG0 

3950-C  ALPH1  IS  A  CORRECTION  FACTOR  FOR  CUMULATIVE  NORMAL  ARCUEMENT 
3960-  ALPH1-1./11.+PALPH) 

3970-  TS-ZL0-'<ZLXALPH1XSIGX> 

3980-  TG-2.XSIGX2XTXTXSICH2'ZL2 

3990-  T1S-ZL02/ZL2 

4000-  T7-T1SXTXTXSICH2 

4010-  GO  TO  £ 

4020-C  SHEAR  DEPENDENT  CALCULATIONS  XXXXXXXXXXXXXXXXXXXXXXX* 

4030-5  IF (SHEARV-OLDSHR)  6,8,6 

4040-6  OLDSHR- SMEARY 

4050-7  T21-SHEARYXSHEARYZ1 .32S97S 

4060-  T8-T6XT21 

4070-  T9-T7XT21/ZL2 

4080-  GO  TO  9 

4090 -C  X  DEPENDENT  CALCULATIONS  XXXXXXXXXXXXIXXXXXXXXXXXXXXX 

4100-8  IF(XX-OLDX)  9,116,9 

4110-9  OLDX-XX 

4120-  X-XXX6080./S28O. 

4130 -C  Tie  INTRODUCES  ASVWCTRY  (MOST  NOTICABLE  FOR  ’O’  UIND  AND  *0*  SHEAR 

4140-C  CONDITION)  UHCN  COUPLED  UITH  TIC  CRITERIA  FOR  Til.  T10  SHOULD 

4150-C  READ!  TIO-ABS(X  )+2.XSIGX 

4160-10  T10-X*2.ZSICX 

4170-  T11-1.  +  (8.XABS(T10)WL 

4180-  IF(Tll-4.)  12,12,11 

4190-11  TIM. 

4200-12  T22-T11XSIC02 

4210-  T30-TSXX 

4220-  IF<T30-6.)  35,36,36 

4230-36  T30-1- 

4240-  GO  TO  37 

4250-35  T38-CUNN0R  ( T30 ) 

4260-37  T12-T9<T10<T19 

4270-43  IFCX)  13,14.13 

4280-14  T13-1. 

4290-  GO  TO  IS 
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■4300-13  1F<ZN-1.»  113,114,113 

4310-114  T 13-EXP (-<A8S(X  >/ZL  >  > 

43a«*  CO  TO  IS 

4330-113  T13-EXP(-(ABS<X)/ZL>*tZN> 

4340-15  SICV2-T22*T8*T12 

43S0-C  STANDARD  CROSSRANCE  DEVIATION 

4360-  SIGY-SIGV2**.S 

4370-C  TIflE  OF  ARR1UAL  OF  FALLOUT  AFTER  BURST 

4380-  TAPR-(.2S+tT15*T16»Tie*T*Tt2.»SICX2)/T14)**.S 

4390-  BETA-ALOC(TARR/31.G> 

4400-  ZLD- “,287- . 52* BET A- .04475* BET At  BETA 

4410-C  COWERS  I  ON  FACTOR  TO  CHANGE  DH+1  HOUR  DOSE  RATE  TO  ABSORBED  DOSE 
4420-  BIO-EXP(ZLD) 

4430-C  G( X )  IS  THE  EXPRESSION  DEFINED  BV  SCHNIDT  ON  PAGE  4 
4440-  GX-l.lSlSlS*IT30*T13)/lZL*T2e) 

44S0-C  THE  FOLLOWING  C(T)  IS  OBTAINED  BV  NULTIPLVING  C(X>  BV  UIND 
4460-  GT-1 . 151515*  (UIND1T1 3  )/'<ZLtT20) 

4470-C  THE  FOLLOUINC  EXPRESSION  IS  GIT)  MULT IPLIED  BV  THE  CUNULATIUE  NORHAL 

4480-C  FUNCTION 

4490-  GT1-GT*T30 

4500-  IF(UIND)  27,27,53 

4510-53  T23-(2.*X)/'(UIND11 .151515) 

4520-  IFCT23-10.1  28,28,27 

4530-27  ALPH22-1. 

4540-  CO  TO  29 

4S50-28  T24-CUHN0R(T23> 

4S60-C  ALPH22  IS  A  CORRECTION  FACTOR  FOR  CROSSRANGE  GAUSSIAN  DISTRIBUTION 
4S70-  ALPH22- 1 . / C 1 . +PALPH*  < 1 . -T24 ) ) 

4580-29  ALPH2-ALPH221ALPH22 

4590-  IF  (  S IGVA2  )  91,90,90 

4S00-C  TO  CALCULATE  V,  GIUEN  X  AND  A  DOSE 

4610-90  S1GVA2-S1GV21ALPH2 

4620-  CO  TO  17 

4630-91  DHXO-T30*T13*T4*FFRAC/SICV 

4640-  DBXO-DKXOSBIO 

4650-  DOSEL- ABS(S1GYA2 ) 

4660-  IF  (DOSED  95,117,95 

4670-95  IF(DHXO/DOSEL-1. )  94,94,92 

4680-94  DH-O. 

4690-  GO  TO  93 

4700-117  DH-0, 

4710-  DB-0. 

4720-  RETURN 

4730-C  THIS  STEP  CALCULATES  VY  GIUEN  AN  INPUT  DOSE  RATE  AND  XX 
4740-92  DH-ALPH22*SIGY*S«RT<2.*ALOG<DHXO/DOSED  >*5280-/6080. 

4750-C  THIS  STEP  CALCULATES  YY  GIUEN  AN  INPUT  DOSE  AND  XX 
4760-93  DI-ALPH22*SIGY*SQRT<2. XABS( ALOG< DBXO/DOSEL ) >  )*S280./60®0. 

4770-  RETURN 

4780-116  IF(SIGYA2>  91.16,16 

4790-C  V  DEPENDENT  CALCULATIONS  **************************** 

4800-16  IF( YV-0LDY1  17,19,17 

4810-17  OLDY-YV 

4820-  Y- YY*6080./5280. 

4830 -C  CROSSRANGE  GAUSSIAN  DISTRIBUTION  FUNCTION 
4840-16  T16-EXP(-.S*Y*V/(ALPH2*SICV2>>/SIGY 

4850-  GO  TO  22 

4860-19  IF  (FFR  AC -OLDER  AC  >  20,22,20 

4870-29  OLDFRAC-FFRAC 

4880-21  DH-DDHtFFRAC 

4890-  DB-DDBXFFRAC 
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4900*  RETURN 

-4910*22  DDH-T3a*T16*T13»T4 

49E0-  DOE-PDHtEIO 

•1930.  00  TO  23 

•49-10*  END 

4550-  FUNCTION  GANNA 4 TN > 

4962-C  C-ANNA  FUNCTION  APPROXlflATED  TRON  HASTINGS  P.1S6 

■4970-  OmHHA-1  .  ■.TMU-6.S766Sa6?*TI1«<e.9?78i?8l*TH»<-0.9235627*Tn»l 

4920-  1 0 . 673990S0*  TH*<  -0. 3262793* TIU0. 07673206  )  )  )  >  > 

4990*  PETURN 

SO0O-  END 

EO10-  FUNCTION  CUNNORITA ) 

S020-C  CUnULATIOE  NORNAL  APPROXIMATED  FRON  HASTINGS  P.  186 
S030-  TM-AfS<TA-l .414213562) 

5040-  7WP-I .  *mt< .  141 13331  *TB*f  .  0866  4027*  T)4t  ( .  02743340* TMTI- .00039446* 

5050-  • TM«. 00328975) )>> 

5060-  TI1P-TMP1X8 

5070-  SUM-l.-l./TUP 

5020-  IFITA.LT.O. )G0  TO  1 

5090-C  THIS  STEP  INTEGRATES  THE  NORMAL  DISTRIBUTION  FUNCTION  FROM 
S100-C  INFINITY  TO  ZERO.  THIS  ALLOUS  A  MAX  I  HUN  AT  GROUND  ZERO. 

5110-  CUNNOR  *.51(1.  *SU11 ) 

5120-  CO  TO  2 

5130-1  CUNNOR 5* <1. -SUN) 

5140-2  RETURN 

5150-  END 
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INITIAL  CONDITIONS! 


CALCULATED  H*1  HOUR  DOSE  RATE  CONTOURS 


CALCULATED  H*1  HOUR  DOSE  RATE  CONTOURS 
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SUBROUTINE  DOSE 


DEFINITION  OF  FORTRAN  TERMS 


Inputs : 

YY 

XX 

SHEARY 

WIND 

FFRAC 

YIELD 


Crosswind  distance  in  nautical  miles 

Downwind  distance  in  nautical  miles 

Crosswind  component  of  shear  in  knots/kilo- 
feet 

Effective  fallout  wind  in  knots 
Fission  fraction 
Yield  in  megatons 


Conversion  from  nautical  to  statute  miles: 


6080 

5280 

= 

1.151515  and  (1.151515)2  = 

1.325975 

Terms : 

YMT 

s 

£n (yield) 

T3 

= 

2  x  10 6  .  yield  +  (R~”x  •) 

(SNC  .  yield) * 

SIGO 

= 

o  -  exp  { .  7  +  ■<-nyeld).  - 
o  3 

3. 25/(4.  + 

( £n  (yield)  +  5 . 4 )  2 ) 

(statute  miles) 

SIG02 

= 

a  2 
o 

(statute  miles) 

Tl 

= 

£n(yield)  +  2.42  (used  in 

"H"  calculation) 

H 

= 

Hc  =  44.  +  6.1  An(yield)  - 

.205  ( in (yield)  +  2.42)  . 

| £n (yield)  +  2 . 42  | 

(kilof eet) 

♦Roentgens  abbreviated  by  "R". 
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SIGH 

SS 

o.  =  .  18H 
h 

SIGH2 

= 

°h2 

T2 

= 

H/60 

T=T 

c 

= 

(1.0573203) (12 (H/60)  -  2 

(1  -  .5  exp- (H/2 5) 2) 

ZLO 

= 

Lq  =  Wind  (T)  (1.151515) 

ZLO2 

= 

L  2 
o 

SIGX2 

= 

o  2 

X 

= 

,  L  +  8o  2 

0  2  (_° _ _o  ) 

o  'L  +  2o 
o  o 

SIGX 

= 

°x 

ZL2 

= 

LZ  =  Lo2  +  2ox2 

ZL 

= 

L 

T14 

= 

L  2  +  .5o  2 

O  X 

,  L  2  +  o  2 
_  _  (  O  X  ) 

n  L  2  +  .5o  2 

O  X 

ZN 

T20 

= 

1. 

or 

T20 

= 

GAMMA(1.  +  1/ZN) 

T4 

2  x  106  (yield) 

L  T(1  +  1)  /2tt 

n 

PALPH 

•001(H) (Wind) (1.151515) 

n 

u 

o 

(kilofeet) 

(kilofeet) 

(dimensionless ) 

(H/60) 2 )  . 

(hours) 

(statute  miles) 

(statute  miles) 2 
(statute  miles) 2 

(statute  miles) 
(statute  miles) 2 
(statute  miles) 
(statute  miles) 2 

(dimensionless) 

(dimentionless) 

■R-statute  miles 
'  hour 

(dimensionless) 
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ALPH1 


=  a  i  = 


1  +  .001(H) (Wind) (1 . 151515) 
a 


L 


(dimensionless) 


T5 

_ 

o 

(per  statute  mile) 

(Laxox) 

2a  2T2a.  2 

T6 

_ 

x  n 

(hr2  kilofeet2) 

L  2 

T15 

— 

o 

L7 

(dimensionless) 

L  2 

T7 

TV 

(hr2  -  kilofeet  ) 2 

T8 

2ax2T2ah2 (Sheary) 2 (1. 325975) 

(statute  miles)2 

T9 

Lo2T2ah2 (Sheary) 2 (1.325975) 

(dimensionless) 

LH 

X 

= 

Downwind  distance  in  statute 

miles  =  x 

T10 

= 

x  +  2a 

X 

(statute  miles) 

Til 

= 

1  +  (8 1 x  +  2ox!)/L 

or 

(dimensionless) 

Til 

= 

4. 

T22 

= 

4a  2 
o 

or 

(statute  miles) 2 

T22 

= 

(1  +  8  |  x  +  2a  | )  /L)  a  2 

X  o 

T30 

= 

Cumnor (1. ) =  $ (1. ) 

or 

(dimensionless) 
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L  x  Lx 

Cumnor  (= — — — )  =  4>(? - — ) 

LajO^  Y  Laiax 


{LoTahsheary(1*151515))  (x  +  2a  >2 

- 57 - 


(statute  miles) 2 


exp  -  (-^-) 


(dimensionless) 


exp  -  (-L^*-)n 


SIGY2 


02=02(1+8IX+  2ox!)  , 

y  O  - - - + 


2(oxTahSheary (1,151515) } 2  + 


( (x  +  20^) LoTa^Sheary (1 . 151515) ) 2 


(statute  miles) 


SIGY 


0y  =  (SIGY 2) 


(statute  miles) 


TARR 


(7-  (X  +  20) 2T2  +  20  2)  . 
‘a  ‘  <'25  +  “ - L  »  ?  .5.  '  1 

O  X 


(hours) 


BETA 


£n(31. 6) 


(dimensionless) 


t  t 

-.287  -  .52iln(3-1^g)  -  .  04475  Unf-jj—)  2) 


(dimensionless) 
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PJWH  m* 


t  t 

BIO  =  exp- (.287  +  .52in(j~)  +  .04475 


T23 


(2x) 


Wind(l. 151515) 


(hours) 

(hours) 


If  T23  =  10  Alph22  =  a2  =  1 


(dimensionless) 


T24 


=  Cumnor  (prr 


2x 


Wind(l. 151515) 


) 


(dimensionless) 


If  T23 
Alph22 


<  10: 


=  a2= 


1  +  (.001(H)  (Wind)  (1. 151515), (1_Cmnnor(  2x  )} 

o  wina 

o 


(dimensionless) 


Alph2 

SIGYA2 


=  a2' 


2  2 
°y  «2 


(dimensionless) 


(statute  miles) 2 


DHXO 


LoX 

cumn°r  (^^-J  exp 


(ill)" 

v  L  ’ 


(2  x  10 6  (yield)  FFRAC)  _  £ 

Lr(l  +  l/n)V2“ay - fx  (R/hr) 


DBXO  =  DHXO  X  BIO 


(R/hr) 


_  ,  V  .  , [x| tn 

Cumn°r  (^7l5_)  SXP  =  (  L  } 


(2  x  10  ) yield (FFRAC) 

LT  (1  +  l/n)/2TT0y  *  810 


DH 


=  H  +  1  Dose  Rate 


(R/hr) 
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DB 


=  Biological  Dose 


(R) 


=  Crosswind  distance  in  statute  miles  =  y 


T16 


e*p->=  (j^r)  =  fy 


(per  statute  miles) 


DDH 


Lx 

Cumnor  (= - — ) 

Lena/ 


exp-%  (  "—^-7  2  ) 

- — <M,n 


(2  x  10 6 ) yield 
LT(l'  +  l/n)/27 


f  .  f 
x  _y_ 


fission  fraction 


H+l 

FFRAC 

(R/hr) 

DDB 

=  DDH  X  BIO 

(R) 

DH 

=  0 

(R/hr) 

or 


Cumnor  ( 


L  x 
o 

Lai  a 


exp-%  (-  X*  z)2  x  10 6  (yield)  FFRAC  exp-(-^J-)n 


a2  0, 


cyLr(l  +  l/n)/2u 


or 


YDH 


Lx  it 

Cumnor  (--  )  exp-%  (-1^-) n  (2  x  10  ) yield  (FFRAC) 

-  y  1,0  l°X  L 

1.151515(2£n(  (Input  DUJ,)  (L)f(l  +  l/n')72Ti  0 

nti  y 


a  2o 


DB 


=  0 
or 


(nautical  miles) 
(ERDS  or  R) 
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DB 


YDB 


=  DDB  .  FFRAC 

or 


«2°y 

1.151515 


(2in< 


DHXO  .  BIO.  . 
Input  Dose 


h 


(nautical  miles) 
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Appendix  C 

Compar i son  of  Terms  Influencing  Crossrange  Dispersion 


This  appendix  contains  graphical  comparisons  of  the 
effects  of  shear  and  torroidal  growth  on  cr^2.  The  torroidal 
growth  and  shear  effects  are  plotted  versus  time.  The  graphs 
contained  within  this  appendix  are  meant  to  supplement  Sec¬ 
tion  III,  Figures  6-8.  As  is  the  case  for  Figures  6-8,  the 
yield  is  10  MT  and  the  fission  fraction  =  1.0.  The  units 

Sthr 1  *  be  abbreviated  as  mph  on  all  graphs  contained  in 

this  section. 


Wind  =  10  mph 

I.OO-3  Shear  =  .3  mph/kilofoot 


r 


S 


»  ot«  g c s3~iim  ainiuis)  /° 


8: 
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Appendix  D 

Derivation  of  Diffusivity  from  Fick's  Law 


Assume  fallout  cloud  has  one  size  particle  with  average 
velocity  V  ( * )  »  Consider  a  differential  volume  in  spheri¬ 
cal  coordinates  dVOL  where: 


N(r,t)  =  particle  number  density  ( 


and: 


Flux  =  F  =  N  .  V  ( 


#particles  -  st.mi. 
hour 


Initial  condition: 


N(r,0)  =  0 


3N(r,t) 

3t 


-Leakage  -  Absorption  +  Source  =  -V 


Jnet 


and: 


j  =  -DV(F)  =  -DV (NV) 
net 


3N  (r ,  t) 
3t 


=  VDV  (NV) 


VDVF  =  V2F 


In  circular  geometry: 


3N(r,t)  _  ^/32N(r,t)  2  3N(r,t)x 

3t  3 F2  r  3r  ’ 

-  N(r,t)  .2  N (r , t) 

-  D3  r  u  3r 


To  solve,  use  Laplace  Transforms.  Let: 

*{N(r,t)  }  h  Sf  (r ,  s) 
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then:  =  sN(r,s)  -  Sf(r,o) 

From  the  initial  condition  N(r,0)  =  0  when: 

.  t  =  0 

s  =  0 


and  let: 

Therefore: 


Dv  =  D  .  V 


.  2D  Sf' 

sSf  =  DVN"  +  — 


where 


ft,  _  dS((r,s) 
N  ~  dr 


or 


Si"  +  -  N'  -  i-Sl  =  0 
r  Dv 


s  assumed  constant  and  the  solution  - 


/-r  /-i 

r,s)  =  gl  e-/Dv  +  e+’V 


C2  is  assumed  0  for  stability.  To  convert  back  to  time, 
form  82  (Ref.  12;  497)  is  used  on  the  following  expression 


\  Ci  ~/d£(/®) 

N (r, s)  =  e  /uv 


where 


Resulting  in:  N(r,t)  = 


r2 

27dv  T  W  e' 
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If  forced  to  resemble  Gaussian  distribution  of  the  form 


1  e-Js(^)2 

7 57c  e  0 


where  r  and  a  have  dimensions  of  length. 


C  = 


St 

2/ttD 


and 


ct2  =  2D  t  (area) 


If  some  initial  radius  present  at  t  =  0, 


of 2  =  oi2  +  ?Dvt  (area) 


Then: 


then: 
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Appendix  E 

Computer  Program  User's  Guide 


The  purpose  of  this  Appendix  is  to  describe/  from  a  user's 
point  of  view,  the  Fortran  computer  code  employing  the  WSEG 
analytical  expressions  to  generate  iso-dose  rate  contours. 
Specifically,  this  discussion  will  concentrate  on  the  Fortran 
code  generated  during  this  thesis  to  reproduce  the  sample  re¬ 
sults  contained  in  Appendix  A. 

The  discussion  is  divided  between  the  main  program  and 
the  subroutine  Dose.  The  complete  code  is  contained  in  Appen¬ 
dix  B  along  with  a  definition  of  terms  for  the  subroutine. 

Sample  output  for  several  options  is  also  included.  Much  of 
what  is  Said  here  has  been  incorporated  within  the  AFIT/WSEG 
program  in  Appendix  B  as  comments  to  aid  the  user  should  this 
thesis  be  unavailable  for  reference. 

Subroutine  Dose 

The  subroutine  Dose  (hereafter  referred  to  as  "Dose") 
contains  in  compact  form  the  coded  WSEG  expressions  discussed 
within  this  thesis.  It  is  nearly  identical  to  the  original 
subroutine  contained  in  Appendix  A  and  obtained  from  Mr.  Ralph 
Mason.  Minor  modifications  to  several  expressions  were  neces¬ 
sary  to  make  the  code  compatible  with  the  ASD  CYBER  74  computer 
system  to  correct  obvious  typographical  mistakes.  User  instruc¬ 
tions  are  contained  within  Dose  as  comments. 

In  general.  Dose  can  perform  two  functions  depending  on 
the  input  parameters.  First,  Dose  can  compute  a  DH+^  or 


Biological  Dose  based  upon  the  inputs  of  effective  wind,  shear, 
fission  fraction,  yield,  crossrange  coordinate,  and  downrange 
coordinate.  This  assumes  a  ground  zero  at  "0"  crossrange  and 
"0"  downrange. 

The  input  parameters  are: 


YY  =  Crossrange  coordinate  in  nautical  miles. 

XX  =  Downwind  (+x)  or  upwind  (-x)  coordinate  in 
nautical  miles. 

SHEARY  =  Average  shear  in  knots/kilof eet . 

WIND  =  Effective  wind  in  knots. 

SIGYA2  =  Any  real  number  greater  than  0.  Not  used 
when  calculating  DH+^  or  Biological  Dose. 


The  call  to  the  subroutine  is  a  standard  Fortran  call  to  a 
subroutine: 


Call  Dose  (DB, DH, SIGYA2 ,YY , XX, SHEARY, WIND, FFRAC, YIELD) 


The  output  parameters  are: 


DB  =  Biological  Dose  in  ERDS 
DH  =  Dh+1  in  Roentgens/hr 


The  second  function  Dose  can  perform  is  to  generate  a 
crossrange  coordinate  based  upon  a  downrange  coordinate  and 
a  given  Biological  Dose  or  DH+^*  Dose  solves  for  the  cross¬ 
range  coordinate  by  solving  the  following  expression  for  Y: 


-Input  Dh+1  =  f, 


exp~hi^) 

/2tt<j 


(33) 
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exp'Js(^H~)2 

or  -Input  Biological  Dose  =  f  .  Bio  .  77no~  y  *34* 

y 

where  f  is  defined  by  Equation  (22) .  Since  the  fallout  pat¬ 
tern  is  nearly  elliptical,  the  crossrange  coordinate  is  either 
(+)  or  (-)  yielding  the  same  result.  In  either  case,  the  in¬ 
put  parameters  are: 

YY  =  Any  real  number.  Not  used  in  this  option. 

XX  =  Downwind  (+x)  or  upwind  (-x)  coordinate  in 
nautical  miles. 

SHEARY  =  Average  shear  in  knots/kilofeet . 

WIND  =  Effective  wind  in  knots. 

SIGYA2  =  -Biological  Dose  or  -DH+^ 

The  call  for  the  subroutine  is  also  a  standard  Fortran  sub¬ 
routine  call: 

Call  Dose  (YDB,YDH,SIGYA2,YY, XX, SHEARY, WIND, FFRAC, YIELD) 

The  output  parameters  are: 

YDB  =  Crossrange  distance  in  nautical  miles  corres¬ 
ponding  to  an  input  Biological  Dose. 

YDH  =  Crossrange  distance  in  nautical  miles  corres¬ 
ponding  to  an  input  DH+^ 

Originally  Dose  did  not  compute  specific  functions  known 


as  g(x),  g(t),  or  <£.g(x)  but  these  expressions  are  used  impli¬ 
citly.  This  program  was  modified  for  this  tlu  is  to  compute 


g{t ),  g(x),  4>.g(t),  and  cumulative  cp.g(x)  and  makes  them  avail¬ 
able  as  output  via  common  statements  where  GtX)  =  g(x),  G(T) 

=  g(t)  and  CUMGX  =  cumulative  <}>.g(x). 


Main  Program 

The  main  program  (WSEG)  contains  two  separate  sections. 
The  first  section  describes  the  purpose  of  the  program,  the 
necessary  inputs  and  the  expected  outputs.  The  second  sec¬ 
tion  contains  three  iterative  loops,  with  appropriate  read 
and  write  statements,  to  do  the  actual  computation. 

This  program  has  been  designed  to  produce  output  such  as 
shown  in  Appendix  A  and  to  provide  g(t),  g(x),  <f> - g  ( t )  ,  or 
cumulative  $.g(x).  In  general,  this  output  is  an  attempt  to 
characterize  the  fallout  pattern  by  describing  several  iso¬ 
dose  rate  contours  in  terms  of  upwind  length,  downwind  length, 
maximum  width,  and  downwind  distance  to  maximum  width.  Also 
included  is  the  maximum  DH+^  at  ground  zero.  The  required 
inputs  in  order  are: 


FFRAC 


IYIELD 


ISHEAR 


IWIND 


Real  number  specifying  fission  fraction 
for  burst. 

Integer  parameter  specifying  the  number 
of  yields  to  be  evaluated. 

Integer  specifying  the  number  of  shear 
conditions . 

Integer  specifying  the  number  of  wind 
conditions . 
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Integer  specifying  output  including  G(T) 
and  time.  If  desired  enter  "1”,  if  not 
enter  M0". 


Integer  requesting  output  of  G(X)  and 
downwind  distance.  If  desired  enter  "1” , 
if  not  enter  "0". 

Integer  requesting  cumulative  G(X)  for 
each  input  dose  rate  condition.  If  de¬ 
sired  enter  "1",  if  not  enter  "0". 

Real  number  specifying  the  downwind  and 
upwind  marching  interval.  The  units  are 
nautical  miles. 

Integer  specifying  which  iteration  the 
write  statements  for  G(X)  and  G (T)  act 
upon.  I.E.,  if  INT  =  10  then  every  tenth 
value  of  G(X)/G(T)  and  distance/time  will 
be  printed. 

Real  number  specifying  the  yield  of  the 
weapon  in  megatons. 

Real  number  specifying  the  effective  wind 
in  knots. 


Real  number  specifying  the  crosswind  shear 
component  in  knots/kilofeet . 


Real  number  specifying  the  D 


H+l 


the  com¬ 


puter  will  use  as  it  generates  the  output 


parameters . 


The  entire  program  including  the  subroutine  Dose  is  in 
the  form  of  a  computer  card  deck.  The  above  input  parameters 
are  read  into  the  computer  via  standard  unformatted  read  state¬ 
ments.  Adhering  to  common  Fortran  procedure,  the  input  para¬ 
meters  are  coded  onto  data  cards  located  immediately  behind 
the  second  multipunch  card  (also  called  an  End  of  Record  Card) 
in  the  computer  deck  which  separates  the  source  program  from 
the  data.  The  information  on  each  card  begins  in  column  one 
as  either  a  real  or  integer  number.  Should  multiple  inputs 
be  placed  on  one  card,  commas  separate  the  individual  para¬ 
meters.  The  following  list  of  data  cards  indicate  the  organ¬ 
ization  of  the  input  data  to  produce  results  such  as  shown  in 
Appendix  B  for  a  one  yield,  wind  and  shear  condition: 


Card 

1: 

FFRAC 

Card 

2: 

IYIELD, ISHEAR, IWIND, IGT, IGX, ICUMGX 

Card 

3: 

XLEN , INT 

Card 

4: 

YIELD 

Card 

5: 

SHE ARY 

Card 

6: 

WIND 

Cards 

7-15: 

DHI* 

Additional  data  cards  will  be  necessary  if  IYIELD, ISHEAR,  or 
IWIND  is  greater  than  one. 

The  output  of  the  program  will  repeat  many  input  para¬ 
meters  along  with  the  specified  output.  The  output  is: 

*Note  that  this  program  is  designed  to  produce  output  such  as 
shown  in  Appendix  B  for  eight  dose  rate  contours. 
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INITIAL 

CONDITIONS 


G(X) 


G(T) 


DWDMAX 

UPMAX 


DBMAX 


RMAXD 

DGZ 

YYMAX 


yield,  effective  wind,  shear,  fission 
fraction  and  step  size. 

Fractional  deposition  rate  of  fallout  per 
linear  mile.  The  units  are  per  nautical 
mile.  Also  included  is  corresponding 
distance  from  ground  zero. 

Fractional  deposition  rate  of  fallout  per 
unit  time.  The  units  are  per  hour.  Ac¬ 
companying  G (T)  is  its  time  coordinate 
in  hours. 

Distance  in  nautical  miles  from  ground 
zero  downwind  to  the  dose  rate  specified 


by  DHI . 

Greatest  upwind  dose  rate  specified  by 
DHI  from  ground  zero.  The  units  are 
nautical  miles.  Note  this  value  may  be 
(-}  or  (+)  depending  on  the  magnitudes 
of  the  effective  wind  and  yield. 


Maximum  DH+1  on  the  hotline  contained 


within  the  total  fallout  pattern  specified 
by  the  minimum  DHI.  The  units  are 
Roentgens/hour . 

Distance  from  ground  zero  in  nautical 
miles  to  DBMAX. 

The  D„  ,  (Roentgens/hour)  at  ground  zero. 
Maximum  crossrange  width  in  nautical  miles 
of  iso-dose  rate  contour  specified  by  DHI. 


l02 


R2MAXW 


=  Downwind  or  upwind  distance  to  YYMAX. 

Units  are  nautical  miles. 

CUMGX  =  Cumulative  G(X)  bounded  by  the  upwind 

and  downwind  range  data  specified  by  DHI. 
CUMGX  calculated  by  trapezoidal  integra¬ 
tion  and  is  dimensionless. 

The  second  section  of  the  main  program  contains  three 
DO-Loops  which  compute  via  subroutine  Dose  those  parameters 
listed  in  the  output.  The  first  DO-Loop  begins  at  ground  zero 
and  marches  downwind  along  the  hotline  calling  Dose  at  each 
location  via  the  call  statement  on  page  (97) .  The  parameters 
DWDMAX,  G (X) ,  or  G(T) ,  RMAXD,  DGZ ,  and  CUMGX  are  computed  for 
each  location.  The  parameter  DWDMAX  is  compared  with  the  value 
determined  from  the  previous  iteration.  The  maximum  DWDMAX 
and  corresponding  RMAXD  are  stored  for  further  comparison  and/ 
or  output.  The  second  DO-Loop  repeats  the  above  process  in 
the  upwind  direction.  Cumulative  <f> . g  (x)  is  computed  by  trape¬ 
zoidal  integration  over  the  pattern.  The  third  DO-Loop  marches 
downwind  along  the  hotline  from  ground  zero  computing  DBMAX  and 
R2MAXW  using  the  subroutine  call  mentioned  on  page  (98)  .  The 
parameter  DBMAX  and  corresponding  R2MAXD  are  compared  between 
iterations.  Maximum  DBMAX  and  R2MAXD  is  stored  for  further 
comparison  and/or  output.  In  all  cases,  the  length  of  the 
iteration  varies  according  to  input  step  size  and  contour 

defining  the  limits  of  the  pattern. 
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As  a  final  note,  it  is  also  necessary  to  preface  the 
WSEG  source  program  with  several  control  cards  in  order  for 
the  computer  to  compile  and  execute  the  program  properly. 

Like  the  data  cards  discussed  earlier,  the  information  is  coded 
beginning  in  column  one  of  each  card.  The  following  control 
cards  represent  the  minimum  required  to  successfully  run  the 
program. 

The  first  control  card  is  the  Job  card  containing  a  three 
letter  identifier,  system  preference,  computer  memory  require¬ 
ment,  computer  access  number,  and  for  AF1T  students,  last 
name  and  box  number.  The  second  control  card  executes  the 
Fortran  compiler.  The  third  control  card  executes  the  binary 
program  generated  by  the  Fortran  compiler.  The  final  control 
card  is  a  raultipunched  End  of  Record  card  which  separates  the 
control  cards  from  the  source  deck.  The  following  is  an  ex¬ 
ample  of  the  control  cards  including  proper  format: 

Card  1:  XXX, STANY, CM60000 . Tllllll , DOE, 0000 

Card  2:  FTN . 

Card  3 :  LGO . 

Card  4:  (7/8/9) 
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(a  )  and  activity  conservation  is  included.  Results  of  the  analy¬ 
sis  of  demonstrate  that  diffusive  growth  is  not  accounted  for  in 
the  model  and  that  crosswind  shear  is  the  dominant,  long  term 
effect.  In  a  comparative  conservation  analysis,  the  WSEG  model 
in  use  today  does  not  conserve  activity  due  to  the  unnormalized 
character  of  the  crossrange  transport  function.  This  effect  is 
substantial  at  yields  less  than  .1  MT.  Activity  not  conserved 
varied  between  31.4%  at  1  KT  and  a  wind  of  60  st.mi.  to  less  than 
1%  at  100  MT  and  winds  of  60  st.mi.  Also  included  is  a  further 
discussion  of  model  limitations  or  inconsistencies  discovered 
either  through  computer  use  during  this  independent  study  or 


during  initial  literature  search. 
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